Conductivity of interacting massless Dirac particles in graphene: collisionless regime 



Vladimir Juricic, 1 ' 2 Oskar Vafek, 3 and Igor F. Herbut 4 

1 Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, 
P.O. Box 9506, 2300 RA Leiden, The Netherlands 
2 Max-Planck-Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany 
National High Magnetic Field Laboratory and Department of Physics, 
Florida State University, Tallahasse, Florida 32306, USA 
^Department of Physics, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6 

We provide detailed calculation of the a.c. conductivity in the case of l/r Coulomb interacting 
massless Dirac particles in graphene in the collisionless limit when uj ^> T. The analysis of the elec- 
tron self-energy, current vertex function and polarization function, which enter into the calculation 
of physical quantities including the a.c. conductivity, is carried out by checking the Ward-Takahashi 
identities associated with the electrical charge conservation and making sure that they are satisfied 
at each step. We adopt a variant of the dimensional regularization of Veltman and 't Hooft by 
taking the spatial dimension D — 2 — e for e > 0. The procedure adopted here yields a result for the 
conductivity correction which, while explicitly preserving charge conservation laws, is nevertheless 
different from the results reported previously in literature. 



I. INTRODUCTION 

The role of Coulomb electron-electron interactions in 
systems described by massless two-dimensional Dirac 
fermions has been a subject of interest for some time^— 
Discovery of graphene, a single atomic layer of sp 2 hy- 
bridized carbon, and more recently of topological insula- 
tors, both of which support such massless Dirac fermions, 
brought this issue into sharp focus. In particular, which 
physically measurable quantities are modified from their 
non-interacting values, and by how much, would al- 
low deeper understanding of the physics governed by 
electron-electron interactions in these systems. 

When weak, the unscreened l/r Coulomb interactions 
are expected to modify the velocity of the Dirac fermions 

2 

as vf — Vp + ^-lnA/fc, where k is the wavenumber 
measured from the Dirac point. This modification of 
the electronic dispersion is expected to lead to logarith- 
mic suppression of the density of states near the Dirac 
point, an effect in principle observable in tunneling exper- 
iments. In addition, the low temperature electronic con- 
tribution to the specific heat should be suppressed from 
T 2 to T 2 / log 2 T, as shown in Ref. H and the strength of 
this suppression is related to the strength of the Coulomb 
interaction. 

The role of Coulomb interaction in a.c. electrical trans- 
port was investigated by Mishchenko in Ref. 0, who orig- 
inally concluded that the a.c. conductivity cf{uj) van- 
ishes as lo — > and the system is a (weak) insulator. 
Were this the case, the interactions would have dra- 
matic effect on transport since the a.c. conductivity 
of the non- interacting system is finite^ i.e., ctq(uj) — 
ne 2 /2h for uj ^> T. This was later argued to be in- 
correct by Sheehy and Schmalian^, and independently 
by the present authors^ using Renormalization Group 
(RG) scaling analysis. While the former presented only 
a scaling argument, without calculating the correction to 
transport, the latter reported on an explicit calculation 



where 

a{w) = aAl + C -^n-^v) (1) 

with the coefficient found to be C = (25 — 67r)/12 ~ 
0.5125. Note that, since e 2 does not renormalize^ii 7 . any 
change in the cutoff in the above expression for the con- 
ductivity may be compensated by a redefinition of the 
Fermi velocity, vf- At small uj the correction vanishes 
and the non-interacting value of a is recovered. At small 
but finite frequencies, the correction scales as 1/| logw|, 
with the interaction independent prefactor determined by 
C. The numerical value of this correction, which can 
be understood as correction to scaling near the Gaus- 
sian fixed point and which is expected to be universal, 
has since been a subject of debate. In subsequent work, 
Mishchenko^ recovered the functional form (pj, which 
gives metallic conductivity at small uj, but argued for 
a different value of C = (19 - 6tt)/12 ~ 0.01254 which 
happens to be much smaller than the one found by us. 
Technically, the difference originated from different reg- 
ularization adopted in the two approaches. The stan- 
dard momentum space cutoff, motivated by the underly- 
ing discrete lattice structure and reported in Ref. was 
questioned in Ref. [ToL where the correction to conduc- 
tivity was calculated using a cutoff on the l/r interac- 
tion, and argued to be the same regardless of whether 
it is calculated using Kubo formula (current-current cor- 
relator) or continuity equation and density-density cor- 
relator. The same calculational procedure was later ad- 
vocated by Sheehy and Schmalian,— who argued that 
unlike hard-cutoff in momentum space, cutoff on the in- 
teraction leads to expressions obeying Ward-Takahashi 
identity. In addition, they claimed the result obtained in 
such way is consistent with the experimentally measured 
optical conductivity, where, surprisingly, no discernible 
correction to the non-interacting value was reported^ 
In quantum field theories, it seems reasonable that if 



two ultraviolet (UV) regularization schemes give differ- 
ent results for physical quantities, then the regulariza- 
tion that is typically chosen is the one which respects 
charge U(l) symmetry, as is the case for chiral anomaly 
in (3+l)-dimcnsional massless quantum electrodynam- 
ics (QED), for instance.— Here we argue that the reg- 
ularization of the electron-electron interaction alone is 
incomplete and cannot serve as a consistent regulariza- 
tion of the theory. We also show by explicit calcula- 
tion that the dimensional regularization used here pre- 
serves the Ward-Takahashi identity, i.e., that it is con- 
sistent with U(l) gauge symmetry of the theory, and 
that, moreover, has the additional advantage of serving as 
an interaction-independent regularization scheme for the 
whole field theory. The interaction correction to the con- 
ductivity within this regularization scheme is calculated 
independently using the current-current and the density- 
density correlators, which both yield the same number 
C = (11 - 3tt)/6 ~ 0.2625 in Eq. (□}, precisely as a con- 
sequence of explicitly preserved U(l) gauge symmetry. 
Furthermore, we show that while the hard-cutoff regular- 
ization in principle violates the Ward-Takahashi identity, 
the original integral expression 9 - for the constant C is in 
fact UV convergent, and when computed with necessary 
care it unambiguously leads to the same value as quoted 
above. 

A comparison with experiment which has been per- 
formed at high frequencies near the cutof f 19 ' 21 (see also 
Ref. I22I ) may be misleading, since the result for the lead- 
ing logarithmic correction to the conductivity in Eq.(p} 
is valid only at frequencies of the order of lmeV, much 
smaller than the cutoff. As the Coulomb coupling con- 
stant in graphene e 2 /vf is believed to be of order one, we 
expect that in this region the interaction corrections to 
different observables, relative to the values in the nonin- 
teracting theory, should be significant. Why the interac- 
tion correction to the conductivity, in particular, appears 
to be small even in the high-frequency regime is unclear 
at the moment. 

Whereas the results in the collisionless limit (lu 3> T) 
discussed here at least in principle follow from a straight- 
forward application of the perturbative renormalization 
group, transport in the collision-dominated regime (u> <C 
T) requires re-summation of an infinite series of Feyn- 
man diagrams. This is easily seen in the non-interacting 
limit where a finite temperature T produces a finite, lin- 
ear in T, "Drude" S— function response in conductivity, 
~ T5(uj). Collisions due to the electron-electron scat- 
tering lead to broadening of the 5— function and clearly 
the result must be non- analytic in e 2 /vf as the interac- 
tion V(r) — > 0. Alternative approach has been advanced 
in Refs. I23II24I where the leading correction is argued to 
be captured by the solution to the quantum Boltzman 
equation with the collision integral calculated perturba- 
tively in the interaction strength. In the clean limit, the 
conductivity in the collision dominated regime is found to 
increase with decreasing T and proportional to ln 2 (T/A). 
Interestingly, experiments on suspended samples at the 



neutrality point^ 5 . find conductivity which decreases with 
decreasing T. Finally, T— linear increase of the d.c. con- 
ductivity observed in small devices 2 ^, has been argued to 
arise from purely ballistic transport^! where conductiv- 
ity grows with the sample size L and temperature T as 
a ~ TL/Hvf- 

The paper is organized as follows: in section II we in- 
troduce the Lagrangian and the response functions, and 
in section III we discuss different regularization schemes 
for massless Dirac fermions. In section IV, we review 
some well-known results regarding the polarization ten- 
sor and the conductivity. In section V, we explicitly 
construct polarization tensor for the non-interacting the- 
ory, and in section VI we consider the same problem for 
the contact interactions to first order in the interaction 
strength and to O(N). We do not discuss the (RPA- 
like) contribution to the order A 2 which, while simple 
to calculate, does not contribute to transport. The main 
results of the paper are presented in section VII where 
we show that the Coulomb correction to the polarization 
tensor is transversal, as well as that the Coulomb ver- 
tex function obeys the Ward-Takahashi identity within 
the dimensional regularization. In this section, we also 
present calculations of the Coulomb correction to the 
a.c. conductivity using both the current-current corre- 
lator (Kubo formula) and the density-density correlator, 
within the dimensional regularization. Section VIII is 
reserved for further discussion of these results and com- 
parison with previous results reported in the literature. 
Various technical details of the calculations are presented 
in the appendices. 



II. HAMILTONIAN, LAGRANGIAN AND THE 
RESPONSE FUNCTIONS 

We start with the Hamiltonian 

H = J d Z3 r-0 t (r)u F cr a p a -0(r) 

+ J d D vd D v'^ ■(r)V(r)F(|r - r'Q^VMr'), (2) 

where we consider N copies of two-component Fermi 
fields ip(r,r) (which therefore has 2N components), the 
momentum operator p a = —ihd a and a a are Pauli ma- 
trices. Operators in the interaction term are assumed 
normal ordered. Hereafter, the Latin letters a, b are 
used only for the spatial indices, while the Greek letters 
fx, v are reserved for the spacetime ones, and summation 
over the repeated indices is assumed. V(r) is the two- 
body interaction potential, which is left unspecified at 
the moment. Later we will consider two different cases: a 
short-range contact interaction V(r) — u5(r) and the 3D 
Coulomb potential V(r) — e 2 /r with e 2 /vp as the dimen- 
sionless Coulomb coupling constant. To simplify the no- 
tation, we will work in the natural units h = c = ItB = 1. 
When the speed of light c does not appear, we will also 



set vp = 1. In our final results we will restore the physi- 
cal units. 

The corresponding imaginary time Lagrangian is 



C — Cq + Ci 



where 



— + v F a ■ p 

OT 



and 



C = J d D r^(r,r 

d^d D r' p(T,r)V(\r - r'\)p(T, r 



(3) 



V(r,r) (4) 



C 



(5) 



where p(r, r) = ip'(r, r) is the density of fermions. 

The quantum partition function can then be written as 
the imaginary time Grassman path integral 28 



Z = J V^Vif) exp \ - drcj 



(6) 



where the inverse temperature factor /3 = l/(fcsT). In 
the sections which follow, the additional imaginary-time 
index on the Fermi field i}){t,v) inside a path integral 
automatically means that they are considered to be co- 
herent state Grassman fields. We will take T — > first 
and then perform the calculations. Note that in light of 
the discussion in the Introduction, taking T — > first 
automatically sets the collisionless limit. 

By the standard spectral representation theorems we 
can first calculate the correlation functions as imaginary 
time ordered products, Fourier transform over time, and 
then analytically continue to find the physical retarded 
(or advanced) response functions*^ Specifically, for some 
bosonic operator O a (t,r) in the (real time) Heisenberg 
representation, the retarded correlation function 

S r a ?(t - t', r, r') = -i6{t - t'){[G a (t, r), 6 b (t', r')] ) (7) 

can be related to the imaginary time ordered correlation 
function 



S ab (T-r',r,r') = -{T T 6 a (r, r)d b (r' , r')), 



where 



d a (r,r) = e^6 a (r)e-^. 



(8) 



(9) 



In the Eqs.j7][8]) the angular brackets denote thermal av 
eraging 

, , 1 



Tr e-P H 



Specifically, the frequency Fourier transforms 

/oo 
dte im S r ^(t,r,r'), 
-OO 

S ab (in n ;r,r') = / dre 40 " T 5 a6 (r, r, r'), 
Jo 



(10) 

(11) 
(12) 



satisfy 

S£f(fi; r, r') = S ab (in n -> n + z0+; r, r'), (13) 

where the bosonic Matsubara frequency is f2 n = 2im/(3 

for n = 0, ±1, ±2, We will use the above relations in 

what follows when we focus on the electrical conductivity, 
in which case the bosonic operator O of interest will be 
either charge density or charge current. 

For completeness we note that for V(r) — the two- 
particle imaginary time Green's function is 

k)V>W, k')) = /3<W (2^) 2 <5(k - k')G k M(14) 

where 



Gk(iw) 



k 2 ' 



(15) 



which will be used extensively in the later sections. 
Strictly speaking, in any solid state system which sup- 
ports massless Dirac particles, the above propagator is 
valid only for wavevectors smaller than some cutoff A, 
which depends on the physical situation. In the case of 
electrons on the honeycomb lattice, the order of magni- 
tude of the cutoff, A -1 , is determined by the requirement 
that the true electronic dispersion does not deviate ap- 
preciably from the conical (Dirac). 



III. REGULARIZATION SCHEMES FOR 
MASSLESS DIRAC FERMIONS 

Since we are interested in the long distance (low fre- 
quency) behavior of physical quantities, we can use the 
above low energy field theory, given by the above La- 
grangian with the corresponding propagators, provided 
that divergent terms in the perturbation theory are prop- 
erly regularized. In the context of high energy physics it 
is also well known that a quantum field theory of Dirac 
fermions needs to be regularized, 20 and typically there 
is no unique way of doing so. Additional requirements, 
usually based on the symmetries of the theory, determine 
what type of regularization should be employed. 

In case of the theory of the Coulomb interacting Dirac 
fermions, we will require that the U(l) gauge symmetry 
must be preserved, or equivalently, that the charge must 
be conserved. As we show below, dimensional regulariza- 
tion introduced by 't Hooft and Veltman 2 ^ is consistent 
with this requirement. Before discussing this regulariza- 
tion scheme, let us briefly review the hard cutoff and the 
Pauli-Villars regularization schemes in the context of the 
fermionic field theory considered here. 



A. Hard cutoff 

The idea of the hard cutoff regularization is to impose 
a cutoff in the upper limit of an otherwise divergent mo- 
mentum integral. Physically, this is due to the k— space 



restriction on the modes which appear in the theory, a 
condition which appears naturally within Wilson formu- 
lation of the RG.— The singular part of the integral then 
appears dependent on the cutoff scale. Although very 
simple to implement, this regularization scheme is known 
to violate C7(l) gauge symmetry of QED, for instance. 20 
Terms that violate the gauge symmetry appear as a 
power of the cutoff scale and must be subtracted in or- 
der to insure that the gauge symmetry is preserved. On 
the other hand, the typical divergent terms appear as the 
logarithm of the cutoff scale. Of course, the cutoff scale 
must not appear explicitly in any observable quantity in 
order for the theory to be physically meaningful. The 
disappearance of the cutoff scale A indeed occurs in the 
calculation of the interaction correction to the a.c. con- 
ductivity within quantum field theory of the Coulomb 
interacting Dirac fermions, as discussed below Eq. ([TJ. 
However, as we show in Appendix D, and as was antic- 
ipated in Reffl8l. the hard-cutoff regularization violates 
the Ward-Takahashi identity. We are thus led to conclude 
that this regularization scheme is in principle not consis- 
tent with U(l) gauge symmetry of the theory. This con- 
clusion notwithstanding, the particular coefficient C from 
the introduction may be written as an integral which is 
unambiguous and perfectly convergent in the upper limit, 
provided the momentum cutoff is taken to infinity after 
all the integrals have been performed (see Appendix H) . 



B. Pauli-Villars regularization 

Another way to regularize divergent self-energy and 
vertex diagrams in QED is to introduce an additional 
artificial "heavy photon". 20 In Euclidean spacetime this 
leads to the following replacement of the photon propa- 
gator 



1 



f 



f 



n 2 + k 2 n 2 + k 2 n 2 + k 2 + m 2 

and the mass parameter M is sent to oo at the end of the 
calculation. Since the additional fictitious particle cou- 
ples minimally to the fermions, the regularization pre- 
serves Ward-Takahashi identities which relate the self- 
energy to the current vertex. However, as such, this reg- 
ularization is unable to render photon polarization dia- 
grams finite. This can be avoided by introducing addi- 
tional Pauli-Villars fermions,— at the expense of making 
the method complicated^ 

In the context of the (2+l)D massless Dirac fermions 
interacting with static (non- retarded) 1/r Coulomb in- 
teraction, the analog of the Pauli-Villars regularization 
is 

II 1 



|k| |k| Vk 2 + M 2 

Physically, this corresponds to cutting-off the short- 
distance divergence of the 1/r interaction, without af- 
fecting its long range tail. This modified interaction pre- 
serves Ward-Takahashi identities relating vertex and the 



self-energy^ but, just as in the case of QED, it fails to 
regularize the polarization function without introducing 
additional Pauli-Villars fermions. Therefore, as such it 
cannot serve as a complete regularization of the theory. 



C. Dimensional regularization 

Originally introduced in the context of relativistic 
quantum field theory, the basic idea of the dimensional 
regularization is to regularize four-momentum integrals 
by lowering the number of spacetime dimensions over 
which the integral is performed. This procedure was in- 
troduced by 't Hooft and Veltman 2 ^ to preserve the sym- 
metries of gauge theories. It also bypasses the necessity 
to introduce Pauli-Villars fermions and bosons. 

Here we employ a variant of the dimensional regular- 
ization scheme in that the frequency integrals are per- 
formed from — oo to +00 while the momentum integrals 
are analytically continued from D = 2 to D = 2 — e di- 
mensions. Such separation of time from space is used 
because in the case considered here the Lorentz invari- 
ance is violated by the interaction terms. A momentum 
integral is therefore calculated for an arbitrary number 
of dimensions D, and expanded in the parameter e. Sin- 
gular parts of the integral then appear as the first-order 
poles in the Laurent expansion over the parameter e, i.e., 
as terms of the form 1/e, and the finite part is the term 
of order e° in this expansion. 

The following I?-dimensional (Euclidean) integrals are 
frequently encountered in this regularization scheme 20 



d D £ 1 



(2ir) D (£ 2 + A)™ (An) D / 2 T(n) A™ - "? 



and 



d D £ 



1 Dr(n-f-l) 



(2tt) d (£ 2 + A)" (4ir) D / 2 2 r(n) 
1 



A" 



(17) 



where T(x) is the Euler gamma function and A > 0. 

Furthermore, Pauli matrices are also embedded in 
D = 2 — e-dimensional space. We thus use the follow- 
ing identity 



cracr^aa = DSoft + (2 - D)a a 5 aiJ 



(18) 



where the sum over the Latin letters a, b, used only for 
the spatial indices, is assumed. The Greek letters /i, v 
are reserved for the spacetime indices. The last term on 
the right-hand side turns out to be crucial for the proof 
of the Ward-Takahashi identity, guaranteed by the U{1) 
charge conservation. This is discussed in later sections. 
In short, the last term in Eq. fTS)) yields the last term in 
Eq. (IA13[) . If the latter were omitted the Ward-Takahashi 
identity would be violated. As elaborated on in the dis- 
cussion section, the same term also accounts for the dis- 
crepancy between the results found in this work, Eqs. 



(|82f . and the result we found previously [Eq. (|G29[) ] for 
the Coulomb interaction correction to the conductivity 
where the last term was omitted. Details of this calcula- 
tion can be found in Appendix E. 



IV. CONSERVATION LAWS, CONDUCTIVITY, 
AND THE STRUCTURE OF THE 
POLARIZATION TENSOR 

In the interest of self-containment, in this section we 
review some well known results regarding response func- 
tions and U(l) conservation laws. Most of these results 
can be found (scattered) in many body - quantum-field 
theory textbooks i 20 ' 28 

In order to calculate the response functions to external 
electro-magnetic fields, it is useful to define the imaginary 
time correlation function 



IV(r,r) = (T T ^(T,r)j„(0,0)), 



(19) 



where the current " three- vector" j'„ is composed of the 
imaginary time density and current as 

j(r,r) = (p(T,r),j(r,r)) 

= (^ t (r,r)-0(r,r),u F -0 t (r,r)o ! -0(r,r)) . (20) 

In this section we temporarily restore to clearly dis- 
tinguish it from the speed of light c used below. 

By fluctuation-dissipation theorem^ the expectation 
value of the electrical current-density operator J(i, r), in 
real time t, is related to the imaginary time correlator 
U^v^n, q). The latter is the Fourier transform (IT21 of 
the tensor defined in Ea. (fl9l) . The expectation value of 
the Fourier transform of the electrical current-density is 
then 



(J Q (fi,q)) = --n a o(iO„^ft + iO,q)$(n,q) 
h 



-U ab (iQ n -> O + i0, q 



A b (fl,q) 



(21) 



The Fourier components of the electric and magnetic 
fields are related to the ones of the scalar and vector 
potentials as 

E a (Cl,q) = i-A a (n,q)-iq a <f>(n,q), (22) 
c 



B(tt, q) = ie ab q a A b (tt,q), 



(23) 



where e ab is completely antisymmetric (Levi-Civita) rank 
two tensor. Using Faraday's law of induction we can 
further relate the Fourier components of the electric and 
magnetic fields as 



e ab q a E b (D,,q) = —B(tt,q 



(24) 



In condensed matter systems with massless Dirac parti- 
cles, propagating with velocity vp , as the relevant low- 
energy degrees of freedom considered here, the (pseudo) 



Lorentz invariance is violated by interactions. If we were 
to consider finite temperature T the (pseudo) Lorentz 
invariance would be violated even in the non-interacting 
limit. Nevertheless, when spatial 0(2) rotational invari- 
ance is preserved, as is the case for problems studied here, 
the general structure of the imaginary time polarization 
tensor is^ 

Il^ifin.q) = n A (^„,|q|)A^ +n B (ift„, |q|)B M „(25) 
where the three-tensors are 

q a qb 



q^ 



B„ 



(26) 
(27) 



The Euclidean three-momenta appearing in the above 
tensors are 

9iiv = diag[-l,l,LV, (28) 
q^ = gij,v(—ittn,<l)v = («O n ,q)^, (29) 
q 2 = q»9p V qu = n 2 n + q 2 - (30) 
The real time continuity equation 



|p + V.J=0 



(31) 



requires that, with our choice of the imaginary time 
"three" -current j = (p,j), the transversality of the 
n Aiy («r2,q) is equivalent to the condition 

(-iQ, q) M LV(if>, q) = LV(zft, q) (-ifl, q)„ - 0. (32) 

Note that this is explicitly satisfied by the expression 
(|25[) . If, in addition, the Lorenz invariance is satisfied, 
II ^ = lis, and there is no need to separate out the spa- 
tially transverse component of the polarization tensor. 



A. Ward-Takahashi identity and vertex functions 



In addition to the condition (I32p . the continuity equa- 
tion (|3"Tjl constrains the form of the vertex function. If 
we define the four-point matrix function 

^(r'-ry-^r-rV-r") = (T tJ> (t, r)^(r', r')VV> r" 

where the imaginary time "three" -current was defined in 
Eq. (|20|l . then we must have 2 ^ 



^(T'-r,T'-v;r-r",r-r") 



(S(t - T")5 D {r - r") - 5(r' - r)S D (r' - r)) 
x5(r'-T",r'-r"). 



(34) 



The above expression relates the exact imaginary time 
four-point function to the exact imaginary time Green's 
function 



G(T,r) = {T T il;(T,r)^(0,0)). 



(35) 



If we rewrite the Fourier transform of 7r M in terms of the 
vertex function A M as 

7r M (k, ioj] k + q, iu + if2) = 

5k(«w)A M (k, «w; k + q, ioj + ifi)^ k+q (iw + zfi)(36) 

then the Ward-Takahashi identity for the vertex function 
A M can be written as 

(— in, q) M A^(k, iuj; k + q, iui + in) = 

(37) 

This identity has to be satisfied order by order in per- 
turbation theory. In what follows, we show that this is 
indeed the case for the interacting theories studied here 
when we adopt the dimensional regularization. 



B. Electrical conductivity 

To relate the polarization tensor to the electrical con- 
ductivity, we simply need to relate the expectation value 
of the current to the electric field. Since we have the 
response to the electromagnetic scalar and vector poten- 
tials, we just need to relate those to the electric and 
magnetic fields. Finally, magnetic field can be related 
to the electric field using Maxwell's equations. As is well 
known, at finite wavevector q and frequency fi, one can 
define the logitudinal and transverse conductivity as the 
proportionality between the induced current and the lon- 
gitudinal or transverse component of the electric field. 

Using Eqs. ([2T]25l l27|) . we find 



(J a (O,q)> = -n A (O + i0,|q|) ., 
n q z 



nq a 



n 



r *(n,q) 



e2 rr in , n I h n2 1aqb A b (Q,q) 

-n4o + *o,lql) q2(q2 _ Q2) e 
^n B( ^ + .o,|q|)(^-M^^ q) 



Furthermore, Eqs. (I22|) - (l24l) imply 

(J a (n ) q)) = ^(0 + iQ ) |q|) q2 ^M^ 



(38) 



E 6 (fi,q) 



(39) 



From the above equations we can read off the longitudinal 
and transverse electrical conductivity 

^M\) = ^-^>n.4(fi + *o,|q|), (40) 



o-±(n,|q|) 



n q 2 - n 2 
e 2 n B (n + io,| q |) 
h in 



(41) 



For q ^ (and ft ^ 0), mi need not be equal a±. How- 
ever, at q = 0, the a.c. conductivities 

o-||(n,q = 0) = er ± (ft,q = 0) (42) 

due to the 0(2) spatial rotational symmetry. 

In the following, we will work solely in the imaginary 
time - Matsubara frequency space, and since we restrict 
ourselves to T — 0, we will drop the subscript n on iVL n . 



V. NON-INTERACTING LIMIT: V(v) = 

For the sake of completeness, and in order to illustrate 
how the general results presented in the previous section 
appear in the specific solvable problem, we first examine 
II M „(zf2,q) in the limit of vanishing V(r). The Fourier 
transform (fT2"]) of the polarization function (IT9l in the 
non-interacting limit is easily shown to be 



H$(«i,q) 

- N f°° — 

, 2ir 



——^Tr[Gk(iu;)(TuGk + a(iu + ifi)ov] 
(Zn ) 

(43) 



where Co is the 2x2 unit matrix. To this end it is useful 
to define the vertex function 



doj 
,2n 



(2n) D 



G , k(iw)cr / xG k +q(iw + ^), (44) 



in terms of which 



n£,>(ifi,q) = -NTr[P„(q,in)a„] 



(45) 



The above expression is divergent at large momenta (UV 
divergent) as is easily seen by counting powers. Note that 
this appears even in the non-interacting theory when cal- 
culating the response functions. As is well known in the 
context of relativistic field theories, this UV divergence is 
unphysical and to obtain the correct answer a regulariza- 
tion is necessary 20 ' 31 The regularization of choice here is 
dimensional regularization which leads to finite expres- 
sions and which is consistent with U(l) gauge symmetry 
of the theory. 

As shown in detail in the Appendix A, using dimen- 
sional regularization, we obtain 



25 M o - 
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(iil + a ■ q)cr M (in + a ■ q) 



n 2 



(46) 



Performing the trace in Eq. (|45j) we find 

n$(ift, q ) = 

N ( -q 2 —idq x _ -inq y 



16^ 2 + q 2 \ -M qy -q x q y q 2 x +n 2 



(47) 



We can write the above matrix more compactly as 

n(°)(ffi,q) = -g^-Ssg^VnM^ (48) 

where we used definitions from Eqs. (|28H3"0")) . The correla- 
tion function (|47ti48|) is explicitly transverse, as it should 
be, and 

(-in,q) M nW(in,q) = nW(in s q) (-*n,q)„ = o. (49) 

From the above equations we find 



4° ) ^,q)=4 0) (zfi,q)--- v ^Tc?. (50) 

Analytically continuing according to Eqs. (|40M1[) . with 
the branch-cut of the -^/i-function lying along negative 
real axis, we find the well-known expression for the 
(Gaussian) a.c. conductivity 

,[°>(n) = ff f(n) = ^. (5i) 

As a side remark, if we were to define TL$} as a cor- 
relation function of a slightly different "three" -current 
(—ip,j), the result obtained directly from Eqs. (|47M8j) 
transforms as tensor under Euclidean 0(3) transforma- 
tions. In real frequencies this is equivalent to relativis- 
tic Lorentz transformations, due to the invariance of the 
non- interacting Lagrangian Cq. 

Therefore, regulating the UV divergences via dimen- 
sional regularization implemented here leads to finite ex- 
pressions which preserve the required U(l) conservation 
laws. The necessary regularization of the "integration 
measure" , as done here via dimensional regularization, 
is independent of the electron-electron interaction V{v), 
as it must be if the non-interacting theory is to lead to 
finite correlation functions. Therefore, as shown already 
by this example, regulating o nly the "momentum trans- 
fer" as advocated in Refs. EMI is clearly insufficient. 



can calculate the interaction corrections to the polariza- 
tion tensor perturbatively in powers of the interaction 
strength u. Such contact interactions certainly consti- 
tute an idealized special case^ Nevertheless, this theory 
has the advantage that one can determine the first correc- 
tion in u to the non-interacting (Gaussian) polarization 

tensor IljS, found in the previous section, explicitly for 
finite q and ft. We can then test the general symme- 
try requirements listed before. The technique of choice 
is again the (variant of the) dimensional regularization 
of Veltman and 't Hooft introduced in section III. Since 
this interaction violates Lorentz invariance we can also 
use this example to study how the difference between 
11^ and TIb arises in such theory. 

It is straightforward to use the Wick's theorem to show 
that in this case, the first order in u, and to O(N), cor- 
rection to the polarization tensor is 

f d D k f°° dcj f d D p f°° duo' 

{Tr[Gk(iw)cr /J G'k+q(ja; + ifl)G p (iu}')a v Gp-q(iuj' — ifl)] 
+ Tr[Gk(«w)cr M Gk+q(«w + iQ)G p (iuj')G] l+Ci (ioj + iQ)a v ] 
+ Tr[G k {iu)a„G k+ q{iLJ + in)a u G k (iLj)G p (iu')}.} (52) 

The last two terms correspond to the self-energy correc- 
tion, while the first one is the vertex correction. Because 
the self-energy for the contact interaction vanishes, 

the last two terms in the Eq. (|52"f vanish as well. The 
remaining term can be written rather succinctly in terms 
of defined previously in Eq. (|44| as 



VI. SHORT RANGE INTERACTIONS: 

V(r) = u5(r) 

While the problem of (2+l)D massless Dirac fermions 
with the contact interactions is not exactly solvable, one 



OT^(*n,q) - uNTr[V„(q,,m)V4-ci,-iCl)].(5i) 

The above expression is manifestly transverse, i.e., it sat- 
isfies Eq. (|3"2")l . as can be seen from Eq. (l4*51) . 



Namely, inserting Eq. (|46j) and performing the traces we find 



OT^(ifi,q)= 



512(^ 2 + q 2 ) 



q 2 (q 2 - 
iflq^q 2 

iVtq y {q 2 



n 2 ) 
-n 2 ) 



iflq x (q 2 - n 2 ) 

q 2 x ( q 2 v ~n 2 ) + (q 2 y + n 2 ) 2 

~q x q y (3Q 2 + q 2 ) 

I 



4 



inq v (q 2 - Q 2 ) 
-q x q y (3Q? 



n 2 q 2 



n 4 



2Q 2 



Finally, the above tensor can be factorized as given by Eqs. (|25H26|) . and we find to first non-trivial order in the 



contact coupling u 

lb 512 

n B (ifi,|q|) = -^V^Tc?+^(0 2 +q 2 ) + O( U 2 ). 
Id 512 

(55) 

Expectedly, the above expression shows that the inter- 
action correction to the polarization functions Ha and 
Hb are different (note the sign difference in front of q 2 ). 
As stated above, the reason for the difference is that the 
contact density-density interaction term u(ip' (r)ip(r)) 2 
breaks the Lorentz invariance of the non-interacting part 
of the Lagrangian. Lorentz transformations in general 
rotate between density and current and we have purpose- 
fully omitted any current-current interaction. 

We can further test the Ward-Takahashi identity (|5T)l 
for the vertex function (|36|) in this example with the short 
range interactions. It can be readily seen that the first 
order in u correction to the vertex vector is 

<5A^(k, ioj; k + q, iui + in) = — ?i'P Al (q, in). (56) 

It follows from the Eq. (gHI) that 

- iQV (q, 1^) + q a 7> a (q, iil) = 0. (57) 

Therefore the Ward-Takahashi identities ([37]) are satis- 
fied, since, as mentioned previously in this section, the 
self-energy correction vanishes to first order in u for the 
short-range interactions. 

Finally, from Eqs. (|4QII4ip . we can infer that the above 
terms correct only the imaginary part of the a.c. conduc- 
tivity, but not the real part. At q = 0, correction is the 
same for the longitudinal and the transverse components, 
and to this order in u we have 

- irei 1 -^)- (58) 

Again, the equality between <7ii(fi) and a±(n) is guaran- 
teed due to the 0(2) rotational invariance of this theory. 
Note also that the fact that the interaction correction 
is proportional to the frequency is implied by the power 
counting at the Gaussian fixed point of the theory, and 
is characteristic for any finite-range interaction^ 

VII. COULOMB INTERACTION: V(r) = e 2 / |r| 

Armed with the above results we now focus on the 
main part of the paper where we study the effects of 
the Coulomb interaction. Unlike in the previous cases, 
we have been unable to find the explicit expression for 
the first order correction to the polarization tensor at fi- 
nite q and fi. Nevertheless, we have been able to show 
explicitly that the first order correction to the polariza- 
tion tensor is transverse, i.e., it satisfies Eq. (|52")) . This 
is shown using dimensional regularization in D = 2 — e 



introduced in section III. Next, we study the first cor- 
rection to the Coulomb vertex function which must also 
satisfy the Ward-Takahashi identity ([57)) . Since in this 
case the first order self-energy is known to diverge loga- 
rithmically, the first order correction to the vertex func- 
tion should also diverge as e — > 0. This can be found 
explicitly in terms of elliptic integrals to order e _1 and 
e° and the identity (|37l) is also explicitly confirmed. Fi- 
nally, we proceed with the calculation of the electrical 
conductivity, first by using the spatial component of the 
polarization tensor at q = but finite n (current-current 
correlation function) , and then by using time component 
of the polarization tensor at finite but small q and finite 
ft. The final results for the conductivity calculated in 
both ways are found to be the same. Specifically, we find 
C = (11- 3tt)/6 in Eq. (J). 

For unscreened 3D Coulomb interactions V"(r) = e 2 /r 
the effect of screening due to dielectric medium is easily 
taken into account by rescaling e 2 in the above formula. 
The 0(e 2 ) and O(N) correction to the polarization func- 
tion is then 

rrTfrf/.r. \ », f d D k f°° duj f d D P f°° du' 

{Vp-k 

x Tr [G'k(iw)(T At G'k+q(ia; + if2)Gp+ q (za/ + in)a v G v {iui')] 

+ Mc-p 

x Tr [Gk(iuj)(T M G k+q (iu; + ^)G p+q (iu/ + iil) 
x G k+q (iw + iQ)a u ] 
+ U k _ p 

x Tr [G k (iw)(T / _ 1 G k+q (iw + ift)ovG k (iw)G p (zcj')]} (59) 
where 

/9wp 2 
d 2 ve^V{v) = fg-. (60) 

Just as in the case of contact interactions, the first term 
in the expression for STlffl corresponds to the vertex cor- 
rection and the last two terms to the self-energy correc- 
tions. Unlike in the case of contact interactions, however, 
the self energy correction does not vanish. The expres- 
sion ([59]) will be used in later sections as a starting point 
in the calculation of the Coulomb interaction correction 
to the a.c. conductivity in the collisionless regime. 



A. Proof of the transversality of SYlj^ within 
dimensional regularization 

Because, as mentioned above, the explicit evaluation 
of ([59"]) at finite q and f2 yields intractable expressions, 
we proceed by first showing that |59|) is transverse, i.e., 
that it satisfies the condition (|32l) , when dimensional reg- 
ularization employed in this paper is used. As such it 
therefore does not lead to any violation of the charge 



conservation, a virtue questioned in Ref. |l8|. By 2D rota- 
tional invariance, this in turn implies that the Coulomb 
polarization tensor can be written in the form (|25j) . 
To prove (|52")l we follow Ref. 18j and use 

-ittcro + q-cr = G k + q (io; + id) - G^ 1 (iu>) (61) 

to find 

d D k d D p f 00 dtudiu' 
{2tt) d (2tt) d J_ O0 2^~2^ 
Vp-k {Tr [Gk(iw)avG k (iw)Gp(ia/)] _ 
Tr [Gk+ q («w + in)(7 I/ G ! k+q(icj + i£l)G p+<1 (iuj' + in)]} 



L)D fry, J 



At this point it is not immediately obvious that we can 
shift the integration variables k and p in the second term 
by q, which if true would readily yield the desired rela- 
tion (I52"j) , since the frequency integral can be shifted. We 
therefore define a function of frequency and two momen- 
tum variables 



£ Pi q(in) 



d D k 



(27T) 



doj' 



D '-co 27T 



V k - p G k+ *(iu'+iQ) (62) 



in terms of which we have unambiguously 
(-ifi,q)„ail#(ifi,q) 

d D k r°° duj 



N 



2tt 



{Tr[G k (^KGkM£ k ,o(0)] 



- Tr [Gk+ q (iw + if2)cr„Gk+q(W + in)£k, q (in)]} • 

To continue, we need to find an explicit expression for 
£ Pi q(if2). Using the identity (|16j) . Feynman parametriza- 
tion 



1 



dyj-T—T-, . a rtHSfl . (63) 



A a B@ T(a)r[P) J *[yA+(l-y)B] a +P' 
for a — (3 — 1/2, and 



dyy<*-\l-yf- 1 

T(a + P) 



(64) 



we find 



E Pi q(«f2) = 



M P + q)r(i-f)r(-±i)iW; 



(4tt)^ |p + q 



2-D 



r(D) 



(65) 

which agrees with Eq.(12) of Ref. [35J. Note that this 
identity shows that within dimensional regularization, 
Ep.q(iri) = £ p+qj o(in). Moreover, in what follows, there 
is no need to shift the integration variable. Rather, since 
the commutator of the self-energy and the Green's func- 
tion vanishes, 



[G k + q (*w + in),£p, q (tn)] = o, 



(66) 



after a straightforward use of the cyclic property of the 
trace and the identity 

/co 
— G k+q (iu; + iO)Gk+ q (ia; + iO) = 0, (67) 
-co 27T 



we prove that 

(-tfi J q)„an$(tfi,q)=0. (68) 

The same procedure as the one used above also leads to 

OT$(tfi,q)(-ifi,q) I/ = 0. (69) 

This proof holds to all orders of e. The regularization 
technique implemented here is therefore perfectly ade- 
quate and does not lead to violation of the charge con- 
servation. 



B. Coulomb vertex and the proof of the 
Ward-Takahashi identity 

Next, we will demonstrate that the dimensional regu- 
larization used here preserves the Ward-Takahashi iden- 
tity for the Coulomb vertex function. This proof is tech- 
nically more involved than the proof in the previous sec- 
tion, but nevertheless, we find it important to present 
its details since our technique is not widely used in the 
community. We show the desired identity to order e _1 
and e°. Most of the technical details are presented in the 
appendices, and in this section we just present the main 
steps of the derivation. 

The Coulomb vertex function to the first order in the 
coupling constant is 

5A M (p, iv\ p + q,iv + ifi) = 7>£(q, p, id) = 

f d D k f°° dw„ _ .. , _ .. 
V (2^p J ^^P-kGkMo>G k +q(zw + ift). 

(70) 

The integrals on the right are logarithmically divergent 
in D = 2 as can be easily seen by powercounting. This 
divergence is related to the divergence of the electron 
self-energy, calculated in the previous section, 



E k («w) = E k ,o(«w) = — cr • k 
o 



7 + ln647r-lnk 2 + 0(e) , (71) 



where the Euler-Mascheroni constant 7 = 0.577 and, as 
before, e = 2 — D. Indeed, if the Ward-Takahashi identity, 

(-in, q)^( q , P , in) = £ p+q (i*/ + in) - e p m, (72) 

is to be satisfied, the vertex function must diverge log- 
arithmically, which manisfests in the dimensional regu- 
larization as the first-order pole in Laurent expansion in 
the parameter e. 

In the second part of the Appendix A we use dimen- 
sional regularization to determine V c (q 7 p, in) to orders 
e _1 and e°. Our final expression for finite q,p and in, 
Eq. (|A18[) . is left as an integral over a Feynman param- 
eter x. We wish to stress that all of the integrals in this 



equation can be performed in the closed form in terms of 
elliptic integrals. However, we found that doing so leads 
to intractable and unrevealing expressions. We therefore 
chose to work with the expression (IA18|) and in effect 
manipulate the integral representation of the elliptic in- 
tegrals. In the limiting case of q = 0, the vertex function 
is determined in the closed form in Appendix A up to, 
and including, e . 

In Appendix B we in turn find that the vertex function 
(|A18|) satisfies 

(-ifi,q)^(q,P,in) = AT(p,q) - e l^wTtf 
x o- (p(r» 2 + q 2 )L(p,q,0)+qM(p,q,r!)). (73) 

Using the dimensional regularization, we then show that 
the function 



N(p,q) = y ( -cr-q+(ln647r-7)cr-q 

- a ■ (p + q) ln(p + q) 2 + a ■ plnp 2 ) 

= E p+q (»i/ + in)-Sp(ti/), (74) 

and, in Appendix C, that L(p, q, SI) = M(p, q, O) = 
0. This proves the Ward-Takahashi identity to the first 
order in perturbation theory. 



C. Calculation of the a.c. conductivity from the 
current-current correlator (Kubo formula) 

In this section we calculate the diagonal spatial com- 
ponent of the Coulomb interaction correction of the po- 
larization tensor, SHxx, at q = and finite ifl. We then 
use this result to calculate the corresponding correction 
to the electrical conductivity. Given the decomposition 
(f2"5)) . one should in principle specify the direction in the 
q-plane along which the limit q — > is taken. For exam- 
ple, if q x is taken to before q v , 5TIxx(i£l, 0) is propor- 
tional to IIb(«^, 0). On the other hand, if q y is taken to 

before q x , then <5n£x(zfi,0) is proportional to II^(zSl,0). 
Similarly, if the limit is taken along a line that forms an 
angle 9 with the q x axis, then SIlixiiQ) is proportional 
to cos 2 0IU(ifi, 0) + sin 2 01Tb (iO,0). However, due to 
the 0(2) rotational invariance, Ha(«^,0) = IlB(jf!,0) 
and the result is independent of 8. We can therefore use 
either Eq. (|40l or Eq. (|41[) along with diagonal spatial 
part of the polarization tensor (j59")l to calculate the a.c. 
conductivity. 

We start by showing that 



5Tl$(q = 0,in = Q) = 0. 



(75) 



This is expected, since a space and time independent 
vector and scalar potential correspond to a pure gauge, 
and as such have no effect on the physics of the problem. 
Within our formalism, this identity can be shown to the 
first order in the Coulomb interaction by first performing 



the integral over the frequencies in Eq. (l59|) . which, as 
can be easily seen, yields 



an<s>(ifi = a >q = o) 



N 



} P 1 



D 



x Tr 



^ C) (0,P,0) 
1 



d D k 



(27T) 



a ■ pavo - • p 



N 



(2ir) D |k| 3 



Tr [(a ■ kcr^ 



p- 

er M er • k) (EkOV 



cr^Ek)] 



Substituting the expression for the self-energy (fTTj) and 
the static vertex (|A20|) . performing the traces and using 
k a kb — > S a bk 2 /D, we conclude that Eq. ([75|) holds, as it 
should. We are therefore free to subtract it from the 
expression for 511$ at either finite O and/or finite q. 

Next, we set q = and consider finite in Eq. (j5"9"1) . 
The polarization tensor can be written as sum of the con- 
tributions from the self-energy correction and the vertex 
correction 

an$(in,o) = anW(*n,o) + fli<2(*fi,o), (76) 

where the self-energy is given by 

n , [°° duoduj' f d D k d D p 

<)(^,0) = 2Nj oo -- J j^nj^V^ 

Tr [Gk(jw)cr Al Gk(iw + i£l)<JuGk(iu>)Gp(iuj')} , (77) 
and the vertex correction reads 

Tr [GkMc^Gk^ + in)G p (iu' + in)a u G P (iuj')$?.8) 

Both of these expressions need to be regulated due to the 
UV divergence. 

In appendix E we calculate both of these contributions 
to the electrical conductivity. The contribution to the 
conductivity coming from the self-energy part expanded 
up to the order e° is found to be 



cr e 



7 — ln(647r) 



The corresponding vertex part is found to be 



<j e 



2 r 



— — 7 + In 647r 



8-3vr 



(79) 



.(80) 



Adding these two terms we obtain the first order correc 
tion to the a.c. conductivity due to the Coulomb inter- 
action 

11-37T 



Sa^ = a a 



(7b 



(i 



-<7 e 



which corresponds to the value 

11-37T 



C = 



6 



(81) 



(82) 



in Eq. ([TJ. We discuss this result in light of previous 
work as well as present day experiments in the concluding 
section. 



D. Calculation of the a.c. conductivity using the 
density-density correlator 

To show that our previous result for the conductivity 
is consistent, we now calculate the longitudinal conduc- 
tivity given by Eq. (j40j) and show that it yields the same 
value of the constant C as in Eq. (1521) . This must be 
the case if the Ward-Takahashi identity and the 0(2) ro- 
tational invariance hold. The longitudinal correction to 

(c) 

the conductivity can be calculated by focusing on oIIqq , 
since -Boo — 0- Unlike for Kubo formula, this component 
of the polarization tensor must be calculated at finite n 
and finite q, since at q = it vanishes. Fortunately, 
we need only the leading order term in the expansion in 
small q 2 , from which we can extract the conductivity. 

According to Eq. (l59l) . the Coulomb interaction cor- 
rection to the density-density correlator reads 

sn$ (in, q ) = (in, q ) + ang? (in, q ) , (83) 

where, just as before, we have separated the self-energy 
contribution 

m («),. n , AT r dudcj' f d 2 k d 2 p 

ankWq) = nJ__ ^ J (27r)2 (2ff)a ^-p 

x {Tr [G k (iu)G p (uj')Gk(w)G k+cl (iLJ + in)} 

+ Tr [G k (iw)G p (w / )G k (w)G k _ q (iw - in)]} , (84) 

and the vertex contribution 

m (b),. n \ »r f°° dujduj' f d 2 k d 2 p 

x Tr [G k (iw)G p (iw')G p _ q (iw' - in)G k _ q (iw - in)] . 

(85) 

The details of the calculations are presented in the Ap- 
pendix F. Here we just state the final result 

( c ) 11 - 37T 2 , . 

oo-Jj = <r e , (86) 

in agreement with the result (|8T|) obtained from the 
current-current correlator. Such agreement is expected 
since, as we have shown to this order in the Coulomb in- 
teraction, the dimensional regularization explicitly pre- 
serves the U(l) gauge symmetry of the theory of the 
Coulomb interacting Dirac fermions. 

VIII. DISCUSSION AND CONNECTION WITH 
PREVIOUS WORK 

Let us now discuss the result ([82]) for the correction 
to the a.c. conductivity due to the long-range Coulomb 
interaction in light of the ones previously reported in the 
literatures^ 

In Ref. HI, the Coulomb correction to the conductivity 
is shown to have the form given by Eq. ([T]), consistent 



with the renormalizability of the quantum field theory 
of the Coulomb interacting Dirac fermions. Moreover, 
the value of the constant C = (25 — 67r)/12 has been 
calculated from the current-current correlator, and us- 
ing hard-cutoff regularization. In Appendix D we show 
that in general hard cutoff violates the Ward-Takahashi 
identity. Nevertheless, it can be shown that the inte- 
gral for C is, despite appearances, in fact UV convergent, 
but sensible to the order of integration. A correct way 
of performing the integral is to integrate both momenta 
up to finite cutoffs, and take the cutoff to infinity after 
all the integrals are done first. This, as shown in the 
Appendix H, corrects the value of the constant precisely 
down to the C = (11 — 3tt)/6. As we showed in Appendix 
G, the previuos result C = (25 — 67r)/12 is also obtained 
when using a version of the dimensional regularization 
in which the Pauli matrices are treated as embedded in 
strictly two spatial dimensions, and which also violates 
the Ward-Takahashi identity, as we argued in Appendix 
D. Technically, the origin of the discrepancy between the 
results for the Coulomb correction to conductivity within 
the two versions of the dimensional regularization may 
be traced if we consider the self-energy correction to the 
conductivity in Eq. ()G2|) and its counterpart with Pauli 
matrices in D = 2 — e, given by Eq. (|E5|) . The dif- 
ference arises from the factor D — 1 = 1 — e which is a 
consequence of the different treatment of Pauli matrices 
within the two schemes. The self-energy piece has a sin- 
gular part proportional to 1/e, and when multiplied by 
a term linear in e coming from D — 1, it gives rise to a 
finite contribution to the self-energy correction. Analo- 
gous situation occurs in the vertex part, and in that case 
the last three terms in the integrand of Eq. (IE13|) account 
for the difference. Namely, when the trace over spatial 
indices of Pauli matrices is taken in D = 2, these three 
terms cancel out, as it may be seen from the term pro- 
portional to (k • p) 2 in the integrand in Eq. (|G6[) . but, 
in fact, when Pauli matrices are embedded in D = 2 — e, 
these terms yield a finite contribution to the conductiv- 
ity, which may be directly checked following the steps in 
Eqs. (|El8)) - (iE23)) . 

On the other hand, in Ref. |lfj, the result for the con- 
stant C — (19 — 67r)/12 has been calculated using three 
different methods, namely, the density-density correlator, 
the current-current correlator and the kinetic equation, 
and it has been argued that in order to obtain the unique 
value for the constant C a short-distance cutoff on the 
long-range Coulomb interaction has to be imposed. This 
regularization is an analogue of the Pauli- Villars regular- 
ization in QED, but without the additional Pauli- Villars 
fermions introduced, that are, in fact, necessary to ren- 
der it consistent^ This value for the constant C has also 
been obtained in Ref. [TH by regulating the short-distance 
behavior of the Coulomb interaction in the same man- 
ner as in Ref. [ljj Although it has been shown that the 
same regularization preserves the Ward-Takahashi iden- 
tity, besides lacking the Pauli- Villars fermions, this reg- 
ularization cannot be applied to the theory of free Dirac 



fermions. Namely, the latter needs to be regularized 
when calculating the polarization bubble. Clearly, this 
cannot be achieved by imposing a short-distance cutoff 
on the long-range Coulomb interaction. Therefore, this 
regularization cannot serve as a consistent regularization 
of the entire field theory. 
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The standard next step when working in dimensional reg- 
ularization is to define new integration variables £ u — 
uj + xfl and £ = k + xq. We then perform the integral 
over £ u] to obtain 



I" 1 f d D £ 

W a) = j o dx J — 

(a-£-xS)a IJ ,(a-£+(l-x)S) 



4VF+A 



A(£ 2 + A) 2 



where we defined 



(A4) 



S = ifl + a ■ q. (A5) 
Since the integration measure is 0(2)-symmetric, only 
the terms even in £ in numerator give a non-trivial con- 
tribution, and we find 



Appendix A: Vertex integrals 

The quantity of interest, which enters into the evalua- 
tion of the bare bubble and the leading order correction 
to the polarization tensor for short range interactions u 
is 

f d^\a f°° dco 
Vn(q,ifl) = J J — G k (iw)^G k+q (ia; + ifl). 

(Al) 



Substituting (|T5|) . using Feynman parametrization ([63] 
for a = /3 = 1, and interchanging the order of integra- 
tions, we obtain 



7^(q, ifl) = dx 
Jo 



d D k f°° dto 



(2tt) d J_ x 2ir [( w + X Q)2 + ( k + xq} 2 + A] 2 

(A2) 



x [iui + a • kjo-p [iuj + ifl + a ■ (k + q)] 



where 



A = x(l - x)(fl 2 + q J ) 



(A3) 



V^(q, ifl) = \ I dx 



(27T) 



D 



VP 



x(l — x)Safj_S 



a M a D(P + A) a 2 (£ 2 + A)i 



(A6) 



We next use the dimensional regularization integrals (|16[) 
and (p~7|) . as well as the identity (fi"8j) to find 



dx 



a u VA-2S fiQ VA 



(ifl + a ■ q)<j fJi (ifl + a ■ q) 



Using Eq. (|A3p . we finally have 



x(l — x) 



A 



V^ifl) 



y/fl 2 + q 2 
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(ifl + a ■ q)^ (ifl + cr • q) 



fl 2 + q 2 



(A7) 



(A8) 



J 



Coulomb vertex 



The Coulomb vertex function to the first order in the coupling constant is 

^ r °° duj 2ire 2 



t,c, [ d k f°° ^ 27re „ t . s „ f 

V^(q,p,ifl) = - j — J ^ | p _ k | Gk(w)q- / »G ! k+q(w + »»)» 



(A9) 



with the free fermion Green's function given by Eq. (|15p . After introducing Feynman parameters using Eq. (I63p . we 
obtain 



V c J<l,p,ifl) 



dx 



d D l 



2ne 2 



(2ir) D \p + xq-£\ 



(a- 1- x(ifl + a ■ q))o>(<7 ■£+(!- x)(ifl + a ■ q)) 



A(l 2 + A)'i 



(A10) 



Now, we consider two terms in the above form of the Coulomb vertex separately. Using the Feynman parametrization 
(|63|) and the D-dimensional integral (|16|) . the first term in the last equation acquires the form 



iD t i 11 rl dy 



(2tt)D \p + xq- t\ VFTA 7T Jo y/y(l-y) J (2*)° I 2 + (1 - v) (y(p + xq) 2 + A) 

(All) 



lT[l-f] f 1 dy 1 



7T (4tt)^ Jo y/y(l-y)-r- (y(p + xq) 2 + A) 1 ^ 
After expanding the integrand to the first order in the parameter e = 2 — D and integrating over y, we have 



(2^)^ |p + xq - ^| VFTA (4tt)t V 2 16 y (p + xq) 2 +x(l-x)(Q 2 + q 2 )y 

(A12) 

The second term in Eq. (|A10|1 . after introducing the Feynman parameter, shifting the variable £ — y(p + xq) — > £, 
and integrating over t, becomes 

d D £ {a- 1- xS) a ll {a ■£+(!- x)S) _ 2 T[l - g] /- 1 rf (1 - yf^ T 



(27r) D |p + xq-^ 2 + A)S 7T (47r)f 7o v 7 ^ (y(p + xq) 2 + A) 1 "* 

2 T[2 - f ] f 1 ^ (1 - y) 5 ^ (ya ■ (p + xq) - arS)^* • (p + xq) + (1 - x)S) 



, dy _ 

7T (4tt)- Jo (y(p + xq) 2 + A) 

2T[2-f] r 1 (l-y)^ ^ 0< r, 



2-# 



„ , dy^ ^ h-, (A13) 

I" (4tt)^ 7o v 7 ^ (y(p + xq) 2 + A) 1 " 



where S — ifl + a ■ q, and we also used the identity (|18j) for the trace of the Pauli matrices over spatial indices. 
Note that the last term in the previous equation arises from the term proportional to e = 2 — D in Eq. (|18[) . If we 
treated Pauli matrices strictly in D = 2, this term would be omitted what would thus lead to the violation of the 
Ward identity, as it may be directly checked in Eq. (|B16[) . Therefore, in order to preserve the gauge invariance of the 
theory, it is crucial to treat Pauli matrices, as well as the momentum integrals, in a general spatial dimension, and 
only at the end of the calculation to take D = 2 — e, and expand in e. After expanding in e, keeping terms up to order 
e , and performing the integral over y in the first term, we have 



I <P£ { «.£- X S M *.£ + (l- X )S) = r^-f] e (pW _ 

J (2n) D |p + xq-^ 2 + A)t (4tt)# "° \^ 2 16 

VA((p + xq) 2 + A) - A e\ 1 



(p + xq) 2 + A 



(p + xq) 2 2 I 4tt 



1 f 1 1 {ya-(p + xq L )-xS)a^ l (ya-(p + xq L ) + (l-x)S) 



2tt 2 Jo y/y(l-y) Vfa + ^q) 2 + A 

r[l-#], 6, (p + xq) 2 / t 



(4tt)t "°r 2'" 16 - Y(p + .Tq) 2 + A 



Sun 1 — — In — e tanh 



{-xSa^a ■ (p + xq) + (1 - x)(p + xq) • ctct^ S)Xi[(p + xq) 2 , A] + a ■ (p + xq^cr ■ (p + xq)I 2 [(p + xq) 2 , A]) 

(A14) 



with A defined in Eq. (IA3|) . The remaining integrals over y read 



X (a,A) 
Xi(o,A) 
X 2 (a,A) 



dy 





y) ya + 


A 


dy 


y 




\/y{i- 


y) ya + 


A 


dy 


y 2 




y/y(i- 


y) ya + 


A 




(A15) 
(A16) 
(A17) 



Therefore, using Eqs. (|A11|) and (|A14|) . we can write the Coulomb vertex in the form 

D 



^(q,P,in) 



A — yj A(A + (p + xq) 2 ) 1 




47T 



as(l -x)m 2 +q 2 ) 



(p + xq) 2 + sc(l - x)(f7 2 + q 2 ) 



( 47r )D/2 y- 2 Hi 

1 — x(l — x)(iQ + a ■ q)a^(iil + a ■ q) 
2^ [ VA((p + iq) 2 +A) 

(— x(zfi + o ■ q)o" M <7 • (p + xq) + (1 - x)(p + xq) ■ aa^iQ + a ■ q)) 



g ■ (p + xq)g M g ■ (p + .xq) 
2(p + xq) 4 



(p + xq) 2 



(p + xq) 2 - 2A 1 - 




(A18) 



where again A = x(l — x)(£! 2 + q 2 ), and we explicitly wrote the functions Xq, X\, and X 2l defined in Eqs. (|A15P - (|A17p . 
The above expression diverges as D — > 2 from below, or equivalently as e — » + . As discussed in the main text, this 
divergence is tied to the divergence of the self-energy and is a consequence of the Ward-Takahashi identity proved 
below. All the integrals over x in the above expression can be performed in the closed form in terms of elliptic 
integrals. (The expressions involving tanh -1 need to be integrated by parts first to bring them to the form easily 
expressible in terms of the elliptic integrals.) However, we found that doing so leads to intractable expressions and 
we thus chose to work with the above form of the Coulomb vertex, in which the integrals over the variable x can be 
thought of as the integral representation of the elliptic integrals. 
At q = the above expressions simplify significantly and we have 



p;(o, P ,ift) 



7 + In 647r — In p z 



21171 



yjQ 2 + 4p 2 



--K 




mi 



-.K 



ym 2 + 4p 2 



(E-K) 



a ■ pa a a ■ p 



vm 2 + 4 P 2 mi 

O 2 (O 4 - 16p 4 )£ + (16p 4 - 2» 2 p 2 )if " 



(A19) 



where the arguments of the complete elliptic integrals of the first and second kind, respectively, are 



K=K 



mi 



V" 2 + 4p 2 



E = E 



mi 



V^ 2 + 4p 2 



Note that to this order, at q = (and at any f2), only the spatial components of are finite. This is a consequence 
of the Ward-Takahashi identity, since to this order the self-energy is frequency-independent. 

Finally, at q = and SI = the integrals in Eq. (|A10[) can be performed for arbitrary D without the necessity of 
expanding in powers of e. The Coulomb vertex then becomes 



This form of the vertex function is used to show that SH^2 (0, 0) = for any D, a fact which is in turn used in the 
calculation of the electrical conductivity. 



Appendix B: Coulomb vertex and the Ward-Takahashi identities in dimensional regularization 

In this appendix we show in detail that to the leading order in the Coulomb interaction coupling constant e 2 and 
to 0(N), the Ward-Takahashi identity, questioned to hold in Ref. 0, is satisfied. As a first step, we define the 
contraction q^V^ = —iVtV^ + loPa- Using Eq. (|A18|) . a straightforward calculation shows that all the terms in the 
contraction proportional to i£l cancel out, and the contraction simplifies to 



? M ^(q,P,^) = / dxl-a-a ) 2J l--n VH ^' -etantr 1 ,/- ^ — H a ■ q 



1 

2^ 



A . Jl-2x)(n 2 +q 2 ) ( A 

- a,f i\h—, — ^mc + a - (p + »q)- — ? x ? 1 - 



(p + xq) 2 +A (p + xq) 2 \ V (p + xq) 2 + A 



A , . , cr • (p + zq o- • qo- • (p + xq 

- — ; -v<7 ■ p + xq)cr ■ qa ■ (p + xq) H ■ r~ 

(p + xq) 4 2(p + xq) 2 



cr • (p + xq)cr • qa ■ (p + xq) / A 3 



(p + xq) 4 V (p + xq) 2 + A 



(Bl) 



In order to show the Ward-Takahashi identity, we first note that the self-energy to the first-order in the Coulomb 
coupling is independent of the frequency. Thus all the terms in the contraction (|B1[) that contain frequency have to 
vanish if the Ward-Takahashi identity holds. In fact, as we will show in what follows, the contraction can be written 
in the form 



e 2 



q^(q, p, ifl) = N(p, q) - - (n 2 + q 2 )W(p, q) - -a- p ^/((l 2 + q 2 ) 3 L(p, q, fi) - -a ■ qv^ 2 + q 2 Af (p, q, ft), (B2) 

where the functions A(p,q), W(p, q), L(p,q, n), and M(p, q, O) are defined in Eqs. ([Bio]) . (|B3|) . (|B17|l . and ([B18]l , 
respectively. This condition is, therefore, satisfied if W(p, q) = 0, M(p, q, £1) = 0, and L(p, q, 51) = 0. Finally, to 
complete the proof of the identity, we will show that N(p, q) = E p+q (w + ico) — £ p (iz/). 
Let us first show that term proportional to ft 2 + q 2 vanishes, i.e., that 

W(p, q)= [ dx ( a • (p + xq) (1 7 2X) - , X ^~% cr • (p + xq)a • qa • (p + xq)) = 0. (B3) 
Jo V (P + zq) 2 (p + xq) 4 J 

Using the identity 

crpcrqcrp = 2pqcrp — p 2 cr • q, (B4) 

we can rewrite the above integral as 

p 2 (l-2x)-(2p-q + q 2 )x 2 



W(p, q) = a ■ p / dx 
Jo 



(P + xq) 4 



cr • q f 1 ( x 2 [2q 2 p ■ q - p 2 q 2 + 4(p ■ q) 2 ] + 2xp 2 [q 2 + 2p q] + p 4 _ 



q 2 Jo V (P + xq) 4 

= ap Wi( p ,q) + f^l W 2 (p,q). (B5) 

q 

When p = q, it is easy to show that both W\ and W2 vanish, and we thus concentrate on the case V = p 2 q 2 — (p • q) 2 > 
0. In order to calculate the integrals W\ and W2, we use the following identities 

K - C dx = p2q2 - 2( ~ P ' q ' )2 - ^ P ■ q ^ q2 q2 C dx CRfi) 
° " Jo (P + xq) 4 2p 2 (p + q) 2 2? 2V J (p + xq) 2 ' 1 > 

Ki= fdx X N4 = q2+P ; 9 q C, dX „, (B7) 







(p + xq) 4 2(p + q) 2 P 2V J (p + xq) 



/ - / 1 A ^ P 2 + P ■ q .PIT dx f R8l 

A2 ^ ' ^^p + xq) 4 ~ 2(p + q) 2 2? + 22?i (p + xq) 2 " 



Straightforward calculation yields Wi(p, q) = and W2(p,q) = 0, and thus W(p, q) = 0. The contraction given by 
Eq. (jBip then simplifies to 



?re 2 e 2 f 1 

9"^(q, P, in) = —a • q Jo - - / dx 



a • q 



2 I -cr • q - a ■ (p + xq) 



a ■ (p + xq)er • qer ■ (p + xq) 

(p + xq) 2 

+ x(l-x)(0 2 +q 2 ) 



(l-2x)(» 2 + q 2 ) 
(p + xq) 2 



a ■ (p + xq)(T • qcr • (p + xq) 

(p + ^q) 4 



sc(l -x)(ft 2 + q 2 ) 



(p + xq) 2 + x(l - x)(ft 2 + q 2 ) 



where 



r[l - f] /- 1 / e (p + xq) 2 _ x f~ 

J ° = ^yWj dx ( v 1 ~2 ln ^ etanh V(F 



x(l-x)(fi 2 +q 2 ) 



+ xq) 2 + x(l - x)(tt 2 + q 2 ) 



After a partial integration in the last term, j7o becomes 
dx ( 1 + e In 4 - e In 



r[i-f] /-i 



v/(p + xq) 2 + x(l - x){Q 2 + q 2 ) + y/x{l - x)(ft 2 + q 2 ) 



(4tt)^ 7 

Finally, after performing another partial integration in the last term of the previous equation, we obtain 



Jo 



r[i-#] f 1 



+ 



{^) D/2 Jo 

e 1 
2 (p + xq) 2 



l-2x 
1-x 



(p + xq) 2 - 2x(xq 2 + p ■ q) 



x(l-x)(0 2 +q 2 ) 



(p + xq) 2 + x(l - x)(0 2 + q 2 



Expansion of the prefactor in e has the form 

which, after keeping the terms up to the order e° in Eq. (|B12I) . yields 



Jo 



2vr 



dx 



1 1, 1 o, x(xq 2 +p-q) 

7 + -[- 7 + ln64.]--ln[(p + q) 2 ] + ^_P_^ 



+ 



2(p + xq) 2 



1 - 2x 



1-x 



(p + xq) 2 - 2x(xq 2 + p q) 



x(l - x)(ft 2 + q 2 ) 



(p + xq) 2 + x(l - x)(0 2 + q 2 



Now, after substituting Eq. (|B14|) into Eq. (|B9|) . the contraction reads 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



q»T c Jq, p, ^) = N(p, q)--(X- p^(Q 2 + q 2 ) 3 L(p, q, n) - -a ■ q-JW + q 2 M(p, q, Q), 



(B15) 



which is, in fact, the form (|B2[) of the contraction, since we have already shown that W(p, q) = 0. The frequency- 
independent part in the above equation reads 



N(p,q) 



=2 /■! 



dx < a ■ q 



2 , i ra i u , \2i 2x(xq 2 +p-q 

7 + In 647r - In (p + qY\ H - 1 

e p + xq 2 



cr ■ (p + xq) a ■ q a ■ (p + xq) 



(p + xq) 5 



while the remaining terms are 
L(p,q,n) = 



1 ^ x 2 (p + q) 2 -(l-x) 2 p 2 
(p + xq) 4 



x(l — x) 



(p + xq) 2 + x(l - x)(il 2 + q 2 ) ' 



(B16) 



(B17) 



and 



M (P) q , o) = i f 1 dx (*\+«r-(i-fp 2 _ 2 \ I *a - «) 



2 7o V (1 -.T)(p + a;q) 2 / V (p + iq) 2 +i(l-i)(fl 2 + q 2 ) 



1 ( x(fl 2 +q 2 )[(3x-2)p 2 + 2x(2a; - l)p • q + x 3 q 2 ] \ / x(l - x) 



(p + xq) 4 J y (p + xq) 2 + sc(l - x)(ft 2 + q 2 ) ' 

(B18) 

The frequency-independent term (|B16|) , after using the identity (|B4j) , and performing the remaining integral, has the 
form 

W(p,q) = y 0ff-q + (-7 + ln64 7 r)(j-q- ( j- (p + q)ln[(p + q) 2 ] + (T-plnp 2 ^ = S p+q (w + iO) -E P (ti/). (B19) 
Here, £ p (izy) is the self-energy defined in Eq. (1711) in the main text. 



Appendix C: Evaluation of the functions L(p, q, fi) and M(p,q, f2) 

In this Appendix we show that the functions L(p, q, f2), Eq. (|B17I) . and M(p, q, Q), Eq. (|B18|) . vanish identically, 
and therefore, the Ward-Takahashi identity is, indeed, satisfied within the dimensional regularization used here. 
In order to evaluate integrals in Eqs. (IB17I) and (IB18|) . we introduce new variables 

p' = P -e-^, w=-, (CI) 

q q 

with p = |p|, q = |q|, and cost/? = p • q/{pq), and thus we can write 

(p + xq) 2 =q 2 (x+p')(x+p'*), (C2) 
with p'* denoting complex conjugate of p' . The integral (|B17|) can now be rewritten in terms of the new variables as 



, _j L q 2 f 1 (l+ P Q(l+p'*)x 3 (l-x)-x(l-x)W* 1 

q *\ L+ wJ ax (x+p'nx+P'*) 2 y/ X (x-i)(x-x + )(x-x-y ( ' 

where x± are the roots of the quadratic equation (p + xq) 2 + x(l — x)(fl 2 + q 2 ) = 0, 

(p + q) 2 -p 2 + n 2 ± V((p + q) 2 -p 2 + n 2 ) 2 + m 2 P 2 

x± = , (C4) 

or in terms of the variables p' and w 



1 + 2Re(p') +w 2 ± J(l + 2Re(p') + w 2 ) 2 + Aw 2 p'p'* . 
X± = 2^ ' (C5) 

Assuming that x+ > 1, we have the folowing sequence 

x+>l>0> x-, (C6) 

which is important for expressing the function L(p, q, ft) in terms of the elliptic integrals, as we will see below. We 
use partial fractions to calculate the integral in Eq. (|C3|) . The first term reads 

x\\-x) = / Vy(i + 2 P ')-(3 + yy»] p' 3 (i+p') 

(x+p') 2 (x+p'*) 2 V S[Im(p')} 3 (x+p') 4[Im(p')} 2 (x+p') 2 '°' 

= -1 + (-^ + — ^ + C.c) , (C7) 
\x + p' [x+p'Y J 



while the second term is 
x(l - xf 



(x +p') 2 (x +P'*) 2 



i(l + p') 2 [p' 2 - 2p'p'* - Re{p')} 
4:{Im(p')} 3 (x +p') 
Ai B 1 



p'(l+p') 



t\3 



A[Im{p')] 2 (x +p') 



+ C.c. 



+ C.c. 



^x + p' (x + p') 2 

We therefore reduced the problem of evaluating the integral (|C3|) to the calculation of the following integrals 



I rn — 



dx- 



1 



1 



(x+p) m ^Jx(x~l)(x-x + ){x-x_y 
with m = 0, 1, 2, and x+ > 1 > > X—. In terms of the integrals I m , the integral L reads 

L = - [1 + 2Re(p')]I + [((1 + p')(l + p'*)A - p'p'*A x )h + ((1 + p')(l + p'*)B - p'p'*B 1 )h + C.c\ 
= -[l + 2Re(p')]I - 



(C8) 



(C9) 



-ip'(l+p')(l + 2p') ip' 2 {l+p') 2 

—h H — rr; — h + <-.c. 



2Im(p') 



2Im(p') 



(CIO) 



The integrals 7 m with m = 0, 1, 2 and a; + > 1 > > x_ have the form (Eqs. (255.00), (255.38), (340.01), and (340.02) 
in Ref. \3Tu 



la 
h 



gF(k), 
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9 

{l+p') 2 a\ 
{2a\k 2 + 2a\ -a\- 3fc 2 )II(a 2 , k))} 



[(a 2 -a 2 )II(a 2 ,fc) +a 2 F(k)] 
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(Gil) 
(C12) 



(C13) 



where F(k), E{k), and II(a 2 , k) are the complete elliptic integrals of the first, second, and the third kind, respectively, 
defined in terms of the corresponding incomplete integrals as Fjk) = F(ir/2,k), E(k) = E(w/2,k), and II(a 2 ,fc) = 
II(7r/2, a 2 , k), with incomplete elliptic integrals defined as in Ref. [33- Here, 



x+ — x_ 

X+(l — X-) ' 
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9 



and 
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y/x+{\ - x~) 
x + +p' 



x+ 1 x+(l+p>) 

Therefore, the integral L has the form 

L = R{p', w)F{k) + Pip', w)E(k) + Im[Gip', w)n(a 2 , k)} 



(C14) 



(C15) 
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Note that imaginary part of a 2 is non- vanishing. In the following, we will show that the coefficient P(p',w) = 0. 
Then, by expressing 7m[G(p', u>)II(p', w)] in terms of the functions Fik) and II(a|, k), with a\ defined below purely 
real, we obtain that the function L vanishes. The coefficient Pip' ,w) reads 
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+ C.c 
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(C17) 



4w 2 (l + Aw 2 ) 
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1 + 2Re{p') - w 2 + v/(l + 2Re{p') + w 2 ) 2 + Ap'p'*w 2 



(C18) 



is purely real, and thus P(p',w) vanishes identically. The coefficient R(p',w) reads 



R( P ',w) = -g(l + 2Re(p'))+g 
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(C20) 



Imaginary part of the product G(p' , ui)II(a 2 , fc) can be expressed in terms of the complete elliptic function of the first 
kind and the third kind, as given by Eq. (419.00) in Ref. 1371 
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where 
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(C22) 



since r 2 + 271 + fc 2 = in our case, with fc 2 given by Eq. (IC14|) . and 01, 61, 71, and 7 2 real. Thus, the integral L 
acquires the form 



where 



and 



L = [R(p' ,w) + Ri(p' ,w)]F(k) + Slip' ,w)n(a 2 ,k), 

g ai{sir 2 - k 2 s 2 m 2 ) + bi(k 2 t 2 m 2 - tir 2 ) 



Ri(p',w) = 



Im(p') 
Si(p',w) = 



m 2 r 2 {s 2 t\ — s\t 2 ) 
9 n 2 (a 1 s 1 - Mi) 



(C23) 
(C24) 

(C25) 



Im(p') s 2 ti - sit 2 

while R(p' , w) is given by Eq. (IC19I) . 

Straightforward calculation shows that a\S\ — b\t\ = and s 2 ti — s\t 2 7^ 0, and thus the coefficient S\[p' ,w) 
vanishes identically. Furthermore, using that a\Si — b\t\ = 0, the form of the coefficient R\{p' , w), given by Eq. (|C24|) . 
can be simplified to 

*^->— 4?5^- (C26> 

Finally, the previous equation together with Eqs. (|C19|) and (|C22I) . since w 2 is purely real, yields R(p' , w)+R\(pt , w) = 
0, and, therefore, the integral L(p,q, Q) given by Eq. (|B17|) vanishes. When the two vectors are (anti)collinear, i.e., 
when Im(p') = 0, the continuity implies L(p, q, 51) = 0. When x + < 1, the analogous calculation shows that 
L(p, q, 57) = 0, as well. 

Let us now turn to the integral M(p, q, O), and consider the case Im(p') 7^ 0, and x+ > 1. In terms of partial 
fractions, this integral reads 

,,2 r 1 r r „'3|-i 1 „/\2 



M 



1 + w z 



dx 



p' 3 (l+P>) 



2ilm(p')(x + p 1 ) 



!\1 



ip' 2 (l+p')(l + 2p') 
2Im(p')(x + p') 



+ C.c. 



1 



\Jx{x - l){x - x + ){x - x-) 2w 
1 

sj x{x — \){x — x+){x — X-) 



- — / dx 



P' 2 {l+P') 
x + p' 



C.c. 



+ p'* 2 + 2(l+p')Re{p')+x~x 2 
2Re\p'{l+p')] + [2Re{p') - l]x + 2x 2 

(C27) 



It follows from Eq. (255.17) in Ref. that 



Ji 

J-2 



xdx 



y/x(x — l)(x — X+)(x — X-) 



[(a 2 -l)IL{a 2 ,k) + F{k)] 



(C28) 



, , * u = 4 i F (k) + 2(a 2 - l)II(a 2 , k) 

o \Jx(x - l)(x - x+)(x - x_) a 



+ 



- 1 



[a 2 £(fc) + (fc 2 - a 2 )F(fc) + (2a 2 fc 2 + 2a 2 - a 4 - 3fc 2 )n(a 2 , fc)] j , 



(C29) 



2(fc 2 - a 2 ) 

which together with Eqs. (IC11[) . (|C12|) and (|C13|) allows us to express the function M in terms of the elliptic integrals 



as 



M = - \A(p' ,w)Il(al,k) + A(p' ,w)*IL(al,k)* +B(p',w)F(k) + X(p' ,w)IL(a 2 ,k) + Y(p' ,w)E(k)] , 
w 

with the coefficients in the above equation of the form 



(C30) 



A( P ',w) =r-.(l-° 



i(l + 2p')(l+w 2 ) 



B(p',w) = 



ip'(l + w 2 ) 
aflm(p') 
1 + w 2 
Im(p') 



2a 

at 



Im(p') 

2 (a 2 - a 2 ) (2a 2 k 2 + 2a\ - a\ - 3k 2 ) 



2(a 2 - l)(k 2 - a 2 ) 
(a 2 - a 2 ) 2 



2{a\ 1) 



+ Re 



2 r? 

a p 



(1 + w 2 )[p'* 2 + 2(1 + p')Re(p')] - Re[p'{l + p')] + 
w 2 (l + a 2 ) 



i(l + 2p'){l + w 2 ) 
Im(p') 



+ w 2 + Re{p') 



X(p',w) 
Y(p', w) 



2a 4 
1 



= 1- 



+ Re(p') + w 2 - 



«; 2 (2a 2 /fc 2 -a 4 + fc 2 -2a 2 ) 



1 + w z 
2Im(p') 



Im 



p' 3 {a 2 -a 2 ) 2 
a 2 {a\ - l)(k 2 - a 2 ) 



2a 2 (k 2 ~a 2 ) 
w 2 (a 2 - 1) 
2a 2 (fc 2 -a 2 )' 



(C31) 

(C32) 
(C33) 

(C34) 



Here, k, a, and a\ are defined by Eqs. (|C14|) and (|C15|) . Since w 2 is purely real, the functions A(p',w), B(p',w), 
X(p',w), and Y(p',w) vanish, and thus the integral M(p, q, ft) given by Eq. (|B18|) is identically equal to zero. This 
result also implies that, because of the continuity, in the case of (anti)collinear vectors p and q, i.e., when Im(p') = 0, 
the function M (p, q, O) also vanishes. When the root x + in Eq, (|C4[) is less than one, the analogous calculation shows 
that the function M = 0, as well. 



Appendix D: Violation of the Ward-Takahashi 
identity within hard cutoff regularization 

In this Appendix we show that Ward-Takahashi iden- 
tity does not hold within the hard cutoff regularization. 
In order to show that, we will follow the same steps as 
in the previous two appendices. Let us first calculate the 
self-energy, given by Eq. (|62|) . 



£ P («w) 



duj' 



d 2 k 2ne 2 Id + a ■ k 



2tt J (2ir) 2 |p - k| uj' 2 + k 2 



ire 



d 2 k 



•k 



(2^) 2 p|p-k| ! 



(Dl) 



which after using the Feynman parametrization (1631) . 
shifting the momentum variable, and performing angular 



integration becomes 



Ep(zw) 



2vr 

• A 

/ dk 



ap dy 



V 



o v 1 - v 
k 



o k 2 + y(l - y)p 2 



(D2) 



where A is the momentum cutoff regulating the ultravi- 
olet divergence of the integral. Remaining integrations 
then yield 



Ep(iw) = — <7-p(lnA-lnp + 21n2) 



(D3) 



Therefore, the divergent part of the momentum integral 
appears as the logarithm of the cutoff which corresponds 
to 1/e pole in the dimensional regularization, see Eq. 



(fTTj) . Note that the divergent parts appear with precisely Analogously, we have 
the same coefficients, but the finite parts are different 
within the two regularizations. 

We now consider the Coulomb vertex function defined 
in Eq. (|70|) which after introducing the Feynman parame- 
ters may be written in the form given by Eq. (|A10[) with 
the integral performed in D = 2. After following the 
same steps as in Appendix A, and performing straight- 
forward integrations, we obtain 



/ 



d 2 £ 1 1 

(2tt) 2 |p + xq - £\ y/p + x{l-x)(n 2 +q 2 ) 



2vr V 2 16 



tanh- J- tf + /2 J • (° 4 ) 

V (p + xq) 2 + x(l -x)(n 2 + q 2 ) J V ' 



d 2 £ {u-l- xS)a IM (a ■£+(!- x)S) _ 6^ f A 1 ^ (p + xq) 2 



x{l - x){n 2 + q 2 ) 



(27T) 



|p + xq-£|(£ 2 + A)t 



2;r 



16 



(p + xq) 2 + x(l - x)(tt 2 + q 2 ) 



VA((p + sq)2 + A )_A\ + 1 { _ + 2 + + 

(p + xq) z I 2tt z k 

(1 - a:)(p + xq) • ct^S] Ii[(p + xq) 2 , A] + <r ■ (p + xcfia^a ■ (p + xq)I 2 [(p + xq) 2 , A] } , 

I 



(D5) 



with Zq, 2i, and 1 2 defined in Eqs. (lAT5|) - (|AT7t . Note 
that the term ~ S^Va is not present in the above equa- 
tion, since the Pauli matrices here are treated strictly in 
D = 2, and thus o a <j^a a = 2<5 M o- Taking the contrac- 
tion q^V^, and following the steps in Appendix B, its 
frequency-independent part reads 



iV(p,q) 



— / dx {a ■ 
4 Jo 



p qj 



fin A + In 4 - lnf 



x(xq 2 



(P 
(P 



xq) 2 
- xq) a 



qer- (p 



2(p + xq) 



_xq)| 5 



■q|] 



(D6) 



which after performing straightforward integrals yields 



A(p,q) = S p+q (ii/ + iO) 



E P (u>) + "g CT ' q ' ( D7 ) 



and therefore the Ward-Takahashi identity is violated 
within the hard-cutoff regularization scheme. The other 
frequency-dependent terms, actually, vanish, since they 
have the same form as within the dimensional regulariza- 
tion, but even if this were not the case they could not can- 
cel purely momentum-dependent term on the right-hand 
side of Eq. (|D7I) that spoils the Ward-Takahashi identity. 
In fact, using exactly the same procedure, one can show 
that within dimensional regularization with Pauli matri- 
ces treated in strictly D = 2 the Ward-Takahashi identity 



is also violated precisely because of the last term on the 
right-hand side in the frequency-independent part of the 
contraction q^T"^. 



Appendix E: Kubo formula and the a.c. 
conductivity within dimensional regularization with 
Pauli matrices in D = 2 — e 



In this Appendix we perform explicit calculation of the 
Coulomb correction to the conductivity within the di- 
mensional regularization, in which both the momentum 
integrals and the Pauli matrices are treated in D = 2 — e, 
which, as we demonstrated, is consistent with the U(l) 
gauge symmetry of the theory, and show that this regu- 
larization yields C = (11 - 37r)/6 in Eq. 

Let us first calculate contribution coming from the self- 
energy part. Using Eq. (|B"5l) and the identity 



<T x a ■ ka x = 2<r x k x - a ■ k, 



the self-energy part of ^Il^^ifi, 0) for /j, 



(El) 



x, 



(SIIosc (ifl, Q), can be written as 



dw 



; r(i-f)rmr(^i) 

(Aiv) D / 2 Y (D) 
k D - 2 



2ir J {2n) D {uo 2 + k 2 ) 2 [{io + Q) 2 + k 2 } 
x Tr{(iu; + a ■ k)[i(a; + f2) + 2a x k x — a ■ k] 
(ico + a ■ k)a ■ k}. 



(E2) 



Performing the trace and the frequency integral, we ob- 
tain 



x r 



' D - 1 



;i-£)r(i-f)rm 



(Att) d / 2 T (d) 



(2ir) D {n 2 +4k 2 ) 2 



(E3) 



After subtracting the zero- frequency part of the Coulomb 
correction to the polarization tensor, and using Eq. (HIT) , 
the self-energy part of the Coulomb interaction correction 
to the conductivity reads 



-32(i Q ne 2 



£)r(i 



2 / 



rmr( 



2 , 



/ dk k 



2D-A 



n 2 + nk 2 

(ft 2 +4fc 2 ) 2 ' 



(E4) 



where er is the Gaussian conductivity of the Dirac 
fermions given by Eq. (|51l) . The remaining integral then 
yields 



a a = — ooe ^ 

2 8-2D (c _ 1)7r 



2 / 



rmr( 



p-r 
2 , 



(4^r(D)r(f; 



(E5) 



cos(Z?7r) 

Taking D — 2 — e and expanding up to the order e°, we 
obtain the self-energy part of the Coulomb contribution 
to the conductivity 



1 



1 



-a e 2 — + - + 7 - ln(647r) + O(e) 



(E6) 



Let us now concentrate on the vertex part of the 
Coulomb correction to the conductivity. Taking the trace 
in 5H x b x given by Eq. (|78l) . we have 



d D P 

(2ir) D 



Vk_p 



[(uj + n) 2 + k 2 }[(uj> + ny 2 + p 2 } 



(w 2 + k 2 )(ui' 2 +p 2 ) 

- [ljoj' + (oj + fi) (u/ + 0)]k • p 

+ 4k • pfc x p x - 2fc 2 p 2 - 2p 2 fc 2 + p 2 fc 2 

- [u(uj' + ft) + lu'(uj + fl)](2k x p x - k • p) 

- u(u + n)(2pg - p 2 ) _ + ft)^ 2 - fc 2 )}(E7) 



Integration over the frequencies then yields 



<mW(in,o) = 2Jv 
i 



d^k 



(2tt) b 7 (2^)° 
[r2 2 (/c :c p a; - k • p) 



fcp(fi 2 +4fc 2 )(f2 2 +4p 2 )' 
4(k-pk xPx + k 2 p 2 -p 2 x k 2 -p 2 k 2 x )}. (E8) 



After subtracting the zero- frequency part of 5H xx (iQ, 0), 
we obtain the vertex part of the Coulomb correction to 
the conductivity 



8a n 
1 



d u k 



D I (o„\D Vk ~P 



(2tt) d J (2tt) d k 3 p 3 (n 2 +Ak 2 ) 

{(k • P k xPx + k 2 P 2 - V \k 2 - P 2 k 2 x )[n 2 



O 2 + 4p 2 
A{k 2 +p 2 )]+Ak 2 p 2 {k-p-k x p x )} 

Cfcl + CT&2 



(E9) 



where 



8cr rj 

1 



O 2 + 4p 2 



d D k r d D P n 2 + A{k 2 + P 2 ) 
J2w)» J J2^)d k - p k 3 p 3 (n 2 + 4k 2 ) 

{(k • pk xPx + k 2 p 2 - p 2 x k 2 - p 2 k 2 x )} , 



(E10) 



and 



Cfc2 = 32<7oil 



d D k 



d D p 



(2tt) d J (2tt) d 
k • p k x p x 

k P (n 2 + Ak 2 )(n 2 + 4 P 2 )' 



14- 



(Ell) 



The contribution Obi, by adding and subtracting Q, 2 in 
the first term in the numerator, may be rewritten as 



(!) , 
fffcl = ^61 + ff, 



(2) , „(3) 



61 



'61 ' 



(E12) 



where 



'bl 



and 



2aM 



4er n Sl 



d D k 


f d D p 


(2n) D J 


(2tt) d 


d D k 


f d D p 


(2n) D J 





k • pk x p x + k 2 p 2 — p 2 x k 2 - p 2 k 


2 

X 




p 

k 3 p 3 











k • pk x p x + k 2 p 2 — p 2 x k 2 — p 2 k 2 



r (2) 



r (3) 
'61 



d D k 



k 3 p 3 
d D p 



(2tt) d J (2tt) 



d D k /• d^p 

WY 5 J 



k 2 



14- 



:§) 2 



A: 2 



V k -p- 



k • pk x p x - 2p 2 fc! 



21,2 



k 3 p 3 



fc 2 + (§) 



(E13) 



(E14) 



(E15) 



The advantage of this decomposition of the integral at is 
that its diverging part is now isolated, and it is contained 
in the integral <J^\\ whereas all the other integrals are 
finite in D = 2. 

We first consider the term aiV . Using the Feynman 
parametrization 



r> + /3 + 7 ] 



dx 



l-x 



A a BPc~< ~ r[a]r[/3]r[ 7 ] Jo 

(1 - x - y) a - 1 x f! - 1 y'>~ 1 
{(1 - x - y)A + xB + yC} a+0+1 1 



dy 



(E16) 



we write 



= — I dx 

7T 



P|k-p|[(f) 2 + fc 2 ] ,, , 

(y-x-yf^x- 1 ' 2 



dy 



(k - xp) 2 + x(l - x)p 2 + y (|) 



3 ■ 



(E17) 



Shifting the momentum k — xp — > k, and retaining terms 
even in k as these are the only non- vanishing ones due 
to the rotational invariance of the integrand, we obtain 



l-x 



= 32a ne 2 I dx I dy (1 - x - y)~ 1/2 x 1/2 



l--)k 2 + x 2 p 2 



to Jo 
d D p p 2 -p 2 x f d D k 



(2ir) D p 3 J (2ir) D 
1 



k 2 + x(l - x)p 2 + y (§)' 



(E18) 



After performing the integral over k, we have 
32er fie 2 / . I 



a W _ 

a w ( 4 7T)^r (f 

(l-x-y) 1 ' 2 dp 



p 



D-2 



{>-i)r(2- 



x(l - x)p 2 + y 



-Vr (3 - #) 



x(l ~x)p 2 + y(%Y 



(E19) 



Integration over p then yields 



(i) _ , uC 2o 2C - 4 2 9 - 2 "r(f-^)r(^i) (i-jy) 



2 ; 

l-x 



bl ' "" (47r)^r( 

x f dxx-^il-x)-^ 1 
Jo 

x {l-x-y) 1 ' 2 . 
Using Eq. and the identity 



dy y 



(E20) 



dy (l-s-yW* = ^^J±(l-x) D -\ 



T(D) 

(E21) 

after integration over y and x, acquires the form 

2o2 D - 4 4 4 - c (i-i>)r(!-£)r(i-f) 



(1) -' 



(47r)^r(f)r(i>) 



(E22) 



Finally, expanding the previous result in the parameter 
e, we obtain 



1 



~ - i (1 + 2 7 - 12 In 2 - 2 In tt) + 0(e) 



(E23) 



Therefore, poles coming from the self-energy and the ver- 
tex parts cancel out, as it should be, since the theory of 
Coulomb interacting Dirac fermions is renormalizable, at 
least to the second order in the Coulomb coupling^ 

Let us turn to the remaining contributions which are 

(2) 

all finite in D = 2. We first consider the term cr^i in 
Eq. (|E14|) . Using the Feynman parametrization (|E16[) . 
we have 



p|k-p|[(§)*+p2] Wo 



dx 



dy 



(p - xk) 2 + x(l - x)k 2 + y (f) 1 



(E24) 



After shifting the momentum p — xk — > p, and integrat- 
ing over p, the term in Eq. (|E14[) becomes 



Ai = -^oe 2 n 3 I dx I dy (l-x-y)- 1 ' 2 



d 2 k 



,-1/2 



c(l-x)k 2 + y(qy fc 2 + (§) 



(E25) 



Integration over the remaining momentum variable then 
yields 



7 (2) 
'bl 



1 



doe 2 dx dy 



o Jo 



1 



y+ \/xy{^- - x) 



yjx(l -x-y) 

(E26) 



After integrating out the Feynman parameter y, the term 



r (2) 



(2) i 2 f 1 i SCC 1 

a,' = -a e / dx 

7T 



o y/x(l - x) 2 a ° ( 7 ) 



We now evaluate the term in Eq. (|E15[) . Using 
Eq. (|E17[) . after shifting the momentum variable, and 
retaining only terms even in p, we have 



7 (3) 
7 bl 



4<r n 3 e 2 / dx dy (1 - x - y) 1/2 x~ 1/2 

Jo Jo 

d D k kl 



{2n)° fc3 

d D p 



k 2 



f i - 2) p 2 - x 2 k 2 



(2ir) D 



p 2 +x(l~x)k 2 +y(%Y 



(E28) 



After integrating over p and setting D = 2, we obtain 

-1 r l-x 



alf = aotfe 2 -^ I dx I dy (1 - x - y) i/2 x 



1 



8n 2 J 



a/2 -1/2 



dk 



fe 2 +(|) :r(l-x)fc 2 + y(§)' 
x 2 fc 2 



2 ar(i - x)fc 2 + y (qy 

Integration over k then gives 



r (3) _ 
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1 rl-x 



cr hl = -— er e / dx dy (1 - x - y) 

4?r Jo Jo 



(E29) 



1/2 



x (l-x)- 1 / 2 V \ (E30) 

^y[jy+Jx^x)} 2 

which, after integrating over the remaining variables, 
yields 



4? = ^ + 3n)a e 2 



1 



(E31) 



Let us now calculate the term avi given by Eq. (|E11I) . 
Using Eq. (|E24|) . shifting the momentum p — xk — > p, 
and retaining only terms even in p, we have 



1 fl-X 

2o / J„ / M. 1 1 ™ „ ,1-1/2 -1/2 



(Tb2 = 4(Toe f2 / da:/ dy (1 — x — y ) ' x 
Jo Jo 

d D k x(k 2 - kl) 



k 

d D p 



k 2 + d) 



(27T) 



p 2 + a(l - x)k 2 + y 



(E32) 



After the integration over p and setting D — 2, we obtain 



C&2 



^di f 1 * dy (l-x-y)- 1 / 2 x 1 ' 2 



k 2 



(2tt) 2 Jo Jo 
x dk 



fc 2 + (§) Z(l-X:)fc 2 + 



(E33) 

Integration over k then gives 

= 2 ^ [(l-^l-x-y)]-* 



Ofc2 



cr e 



cr e 



. dx dy 

47r 7o 7o y/x(l -X) 



1 
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1-z 



■ sec 



(E34) 



where sec 1 a: is the inverse function of sec a = 1/cosrr. 
Finally, integration over x yields 



4 - 7T 



Cb2 



-a e 



(E35) 



Therefore, using Eqs. flES)). (jE23| . (jE27l) . (jE3Tj) . and 
(|E35|) we obtain the first order correction to the a.c. con- 
ductivity due to the Coulomb interaction 

ftr W =a a + + aif + off + CTfc2 = ±^a e 2 , 

(E36) 

which corresponds to the value 

11 (E37) 



C 



Eq. Q. 



Appendix F: Longitudinal conductivity using 
density-density correlator 

In order to obtain the longitudinal conductivity, we 
expand #IIqq (in, q) in Eq. (|53")) , which is the time com- 
ponent (/i = z/ = 0) of the Coulomb correction to the 
polarization tensor (|5§|) . to the order q 2 . 

Let us first consider the self-energy part (fMf which 
may be written as 

4' q) = ^ (in, q ) + sn&v (-in, - q ), (fi) 

where 

xttM^-o \ iu /" rf2k f duJ f d2 P f duj ' 
Vk_pTr [Gk(iw)Gp(«cj')Gk(«w)Gk+q(iw + id)] . 

(F2) 

Using Eq. (|65j) to integrate over the momentum p, and 
taking the trace over Pauli matrices, we obtain 



du> 



d D k 



[Air) D / 2 T (D) 
k D - 2 



2ir J (2ir) D {uj 2 + k 2 ) 2 [(uj + SI) 2 + (k + q) 2 ] 
x [(-uj 2 + k 2 )k- (k + q) - 2uj(uj + n)k 2 }. 



(F3) 



Note that above expression contains a part divergent in 
D = 2 that arises from the self-energy (pj|) . and when 
multiplied with the remaining terms of the order e gives 
a finite result, i.e., the final result does not have a pole 
in e . Therefore, in order not to overlook this subtle 
cancellation, we have to perform the integrations in D- 
dimensions first, and only at the end of the calculation to 
take D — 2 — e, with e — > 0. Expanding the q-dependent 
term in the denominator to the quadratic order in q, and 
keeping only the terms quadratic in q in the expression 



for (JIIqq 1 '' (ift, q), we find 



r(i-f)rmrW) 
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[H)i(a¥] 
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and 



5ri a 12} ( i n, q ) = 2Ne 2 - {l ^) r (^) r v > 
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4(k • q) 2 
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(F6) 



We first consider the term OTqq 11 ' in Eq. (|F5|) . After 
integrating over u> and using the rotational symmetry of 
the integrand, we have 



suil n) (in,c0 



d D k 

D J (2tt) d> 



„D-3 



Av2 r(i-f)rmr(V) 

(Air) D / 2 T (D) 

32/t 4 - \2k 2 n 2 - n 4 



{n 2 



■„2\3 



(F7) 



After integrating over k, and expanding the result in e, 
we obtain to the order e° 

^"("■^-a^-i^M- <F8) 

The term 5H(^ 2 ' (in, q) given by Eq. (|F6I) . after integra- 
tion over the frequency and using the rotational symme- 
try of the integrand, acquires the form 

6U^ 2 \in, q) = -Ae 2 - (1 ~ ^ r ^ r 



(47r) D / 2 r (D) 



d D k 



k 



D-3 



2D J (2tt) d (n 2 +Ak 2 f 
x [16(D - 5)fc 4 + 2Ak 2 n 2 + (3 - D)n 4 } 



(F9) 



Integration over k in the last expression, and expansion 
in e then yield 

mSfw*— i^jg,— (fio) 



Therefore, using Eq. (IF4[) . we have 

5 ^ m ^=k Ne %rr° e *w (F11) 

which together with Eq. (|F1[) gives for the self-energy 
part 

<m^(zfi,q) = ia oe 2 ^. (F12) 
I 



Let us now concentrate on the vertex part of the 
density-density correlator given by Eq. (|85l) . In order to 
calculate this contribution, we need the following integral 
over the frequency 



dui' 

— Gp(iw')G p _ q (iw' - ifl) 

ZTT 



2 (|p-q|+pr + ^ 2 



Hp — q| — p + itta 



2p q 

|p-q| 



+ (—iQa ■ p + p 2 - p • q iosp x q) 



Ip - q| +p 
|p - q|p 



(F13) 



Since the frequency integrals in Eq. ((85J) are decoupled, we use the above equation to perform them separately, and 
then take the trace using the identity 



Tr[(a + b ■ a)(c + d • a)] = 2ac + 2b ■ d. 
Integrations over u> and ui' in the vertex part, given by Eq. (pp)) . yield the coefficients 

i\ i, \f, p • (p - q) A /ii _i , k -( k -q) 

a = -(|p - q| + p) 1 ; I , c - '• -' 



(F14) 



pp-q 



-Qk-ql+A;) 1 



p q p\ ., s p-q+p 

j 1 ."HP x q) z - : — 

p-q v p-qp 



, d 



in 



*|k-q| 

k q k\ |k-q| + / 



(F15) 



in Eq. (|F14[) . where (p x q) 2 = p x q y — p y q x . Notice that 
b and d are three-dimensional vectors. Expanding the 
expressions in Eq. (|F15[) to the order q 2 , we obtain 



q 2 p 2 - (q ■ p) 2 
p 6 

q 2 k 2 - (q-k) 2 
P ' 
P(P • q) - qp 2 



b 
d 



if] 



' 2z (pXq) - 



P 



■ n k(k.q)-q^ .(kxq) 



A- 3 



(F16) 



thus ac = C(q 4 ), and therefore does not contribute to 
the conductivity, whereas 



bd 



2 LP(P • q) - qP 2 ] • [k(k • q) - qfc 2 ] 



p 3 k 3 



(P x q) z (k x q) z 
pk 



(F17) 



Setting D = 2 in the momentum integrals, since there 
are no divergent subintegrals in the vertex part as it was 
the case in the self-energy term, the vertex part of the 



density-density correlator to the order q 2 has the form 



<m$ (in, q ) = sn$? (in, q ) + 5n^ 2) (in, q ), 



where 



ctt(h)/.« n *^ 2 /" rf2 k /" d 2 p 



(p(p • q) - qp 2 ) • ( 


k(k-q) -qfc 2 ) 


p 3 /c 3 


(2 P ) 2 + n 2 




(2fc) 2 + n 2 





d 2 k (k(k-q)-qfc 2 ) •Ii(k,q,n) 



(27T) 



fc 3 ((2/c) 2 + n 2 ) 



(F18) 



and 



p • kg 2 — p q q k 



iV 



(2fc) 2 + n 2 [(2p) 2 + n 2 ] 

ri 2 k J 2 (k,q,n) 

(2tt) 2 fc[(2fc) 2 +n 2 ]' 



(F19) 



Here, we defined 
IxCk.q.n).-? / (F20) 



8 7 (2tt) 



p 2 +(§: 



and 



l a (k,q,n) = - 



1 /" d 2 p p • kg 2 — p ■ q q ■ k 



2 7 (2?r) 



P 



p 2 + (§; 



(F21) 

We consider the term 511^ (zf2, q), given by Eq. (|F18[) . 
and, as a first step, calculate the integral Ii(k). Using 
the Feynman parametrization (|E16|) , this integral can be 
written as 



We use this result to calculate integral over k in Eq. 
(|F19|) . After performing straightforward integrations, we 
find 



Ne 2 q 2 '•' '■' ' 
64^TT| J , 7|1 



SU^^cO^irjz^ I dx I dy 



y/l-x-y 



x(l -x) + y/x(l - x)y 64|S7| 
e 2 q 2 4 - 7r 



(F26) 



Using Eqs. (IF241 and (|F26]l . the vertex part of <5IT 00 
reads 



Il(k,q,fi) 
d 2 p 



-e 2 Q 2 dx dy 
Jo Jo 

p(p • q) - qp 2 



V 1 - x - y 



(2?r) 2 [(p - x k) 2 + x(l - x)k 2 + |0 2 ] : 



(F22) 



which after shifting the momentum variable p — xk — > p, 
and integrating over p, acquires the form 

l 1 (k,q,0)=-e 2 ^ t dx f 1 '* dy ^ X ~ X -y- 



8tt 



a; 2 k(k • q) — qx 2 k 2 
x(l - x)k 2 + |02 + ( X {1 - x)k 2 + fft 2 ) 2 

(F23) 

Using the previous result for the integral Ii (k, q, fi), after 
integration over k, we obtain 



^Ve 2 q 2 f 1 f 1 -* y/l-x-y 



o Jo 



xy 



x 2 



-TV- 



e 2 g 2 



64tt|S1 
e 2 g 2 1/1 7T 



47T 

h 7T I 7T 

3 V 3 



Ifil 16 V 3 4 



1 7T 
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W^VS 4 



(F24) 



We now evaluate the term <$II 00 '(iQ, q) in Eq. (|F19I) . 
First, we compute the integral ^(k, q, Q) in Eq. (|F21[) 
using the Feynman parametrization (|E16|) 



J 2 (k,q,fi) = e 2 f dx f 
Jo Jo 



dy 



1 



y/xy/l -x-y 



d 2 p 



p • kg 2 p • qq • k 



(2tt) 2 ((p - xk) 2 + x(l - x)k 2 + fft 2 



if ft* 

8t Jo Jo 



dy 



) 2 
A- 2 



y/l-x-y [x(l - x)k 2 + f ft 2 ] ' 

(F25) 



OT ooWq) = -jnj — g— . 



(F27) 



which together with the self-energy part (|F12|) yields up 
to the order g 2 



m W,. n v eVH-37T 



(F28) 



Finally, using Eq. (j40|) . we obtain the Coulomb interac- 
tion correction to the longitudinal conductivity 



11 - 3tt 



G 



-<7 e 



(F29) 



in agreement with the result (|8T|) obtained from the 
current-current correlator, which is expected, since the 
dimensional regularization explicitly preserves the U(l) 
gauge symmetry of the theory of the Coulomb interacting 
Dirac fermions. 



Appendix G: Kubo formula and the a.c. 
conductivity within dimensional regularization with 
Pauli matrices in D = 2 



In this Appendix we present the calculation of the 
Coulomb correction to the a.c. conductivity using di- 
mensional regularization, but treating the Pauli matrices 
in D = 2 spatial dimensions strictly . As we commented 
earlier, this leads to the violation of the Ward identity, 
and thus it is incompatible with the U(l) gauge sym- 
metry of the theory, but yields the number obtained in 
Ref. We use Eq. (14"TT) with the Coulomb correction 
to the polarization tensor given by Eq. (I75|) . Since the 
system of Dirac fermions interacting only via the long- 
range Coulomb interaction is isotropic, translationally 
and time-reversal invariant, the trace over spatial indices 
of the polarization tensor at zero momentum and a finite 
frequency is n aa (zf2, 0) = DUb{i^,0), which we use to 
calculate the Coulomb correction to the conductivity. 

Let us first consider the self-energy part obtained from 
Eq. (ITT)) . Using that in D = 2, OaU^Ua — 25 ^, taking 
the trace over the Pauli matrices, integrating over the 



frequencies, and subtracting the zero-frequency part, the 
self-energy contribution reads 



with a a given by Eq. (|E5I) . This result, after expanding 
in e, reads 



8 f d D k f d D p 2ire 2 
a °--D a ° n J (2^Fi (2^p |k-p| 

k- P n 2 + uk 2 
x ~kfk 2 (n 2 + 4fc 2 ) 2 ' 

Performing the momentum integrals, we have 

08— 4_DtV -D+1 U1 M 

a a = —doe si 



(Gl) 



1 1 



~°oe 2 - + -+7-ln647r + C(e) 



(G3) 



The vertex part of the Coulomb correction to the con- 
ductivity is obtained from Eq. (|78l) . The trace over spa- 

n D - 1 r(D/2)T(D)Cos(nD) tial indices of 5n$ in Eq. (|78p at the momentum q = 0, 

using the standard anticommutation and trace relations 
(G2) for the Pauli matrices in D = 2, assumes the form 



, M/ . f°° du f d D k f°° dv f d D p 



(2n)D J_ oo 2n J (2tt) 
((cj + f])(t/ + 0)-2k- gX-a^ + k ■ p) - awk ■ p + fc 2 p 2 
X k ~ p (w 2 + fc 2 )((w + O) 2 + k 2 ){v 2 + p 2 )((v + fl) 2 + p 2 ) ' ( ' 

Performing the integrals over the frequencies u> and v in the above equation, we have 
AnWr-o on ( d ° k I d ° p v k-p(4k-p-» 2 ) 

After subtracting the zero-frequency part of <5II aa (if2, 0), we obtain the vertex part of the Coulomb correction to the 
conductivity 



, d D k f d D p 2Tre 2 k • p [k • p(n 2 + 4k 2 + 4p 2 ) + 4k 2 p 2 ] 

° h D a ° n J (27p J (27F ^-pl k ^(n 2 + 4k 2 )(n 2 + 4p 2 ) = abl + am + f763 ' !(,(,! 



where 



1 o f d ° k f d ° p 2we2 n 2 + 4(k 2 +p 2 ) 

ahX ~ D a ° J (2tt) d J (2tt) d |k-p| kp(n 2 + 4k 2 )(fl 2 + 4p 2 )'- { ' 

diverges in D = 2, as we will show in the following, and the remaining integrals 

. = 32 f^^[^p_2^ kp 

and 

- =1 Q [ d ° k [ d °P 2ire2 [2(k- P ) 2 -fcV][» 2 + 8fc 2 ] 

°b3 D <?o J ^n) D J ( 2?r )D | k _ p | k 3 p 3 {Q 2 +4k 2 )(n 2 +4p 2 ) 1 ' 

are finite in two dimensions. The above decomposition is obtained from Eq. (|G6[) by adding and subtracting J7 2 in 
the term in the numerator multiplying (k • p) 2 . We now further decompose the integral a^i in order to isolate its 
diverging part 

*bi=4i+4?- (G10) 



Here, the term 



_(i)__8 f d D k f d D p 2ne 2 1 

°bi D a & 1 J (2n) D J (2n) D |k - p| kp(Q 2 + 4k 2 ) ' 1 ' 



has the pole in Laurent expansion in e, while the remaining one 



5g } = ~^ J 



d u k f d u p 2ne 



1 



{2ix) D J (2ir) D \k-p\kp(n 2 + 4k 2 )(n 2 +4p 2 ) 



(G12) 



is finite in D = 2. 

We first consider the term divergent in D = 2, namely, cr^. Integration over p in Eq. (jGllI) . after performing the 
standard steps, yields 



d D P i _fc c - 2 r(i-#)[r(^i)r 



(2ir) D p\k - p| 7T (47r)^/ 2 r(Z3 - 1) ' 
while after integrating over k, we find 



The previous expression, after expanding in e, reads 



2 4D - 7 7r D DT(D-l)T(^ 



1 



'51 



1 1 

7+2 



7 + ln(647r) + O(e) 



Thus the poles in a a , given by Eq. (|G3|) . and cancel out, and 



(G13) 



(G14) 



(G15) 



(G16) 



The remaining integrals are finite in D — 2, and can be calculated as follows. Consider the term a b 2 given by Eq. 
G8I) . Using Feynman parametrization (IE16|) . we write 



1 



1 



= — dx dy 



o Jo 



l-x 



{l-x-y)- l ' 2 x- l l 2 



(p - xk) 2 + x(l ~ x)k 2 +y{%Y 



p\k-p\[{%) 2 + P 2 } 7T 

As usual, we now shift the momentum p — xk — > p, and integrate over p, to obtain 

-1/2^1/2 

JO JO 

d 2 k 



(G17) 



I r 1 /-l— a; 

&b2 = —a Q e 2 fl / dx dy (1 — x — yY^^x 1 

27T Jo Jo 



(2n) 2 



x(l-x)k 2 +y(%y k 2 + (%) 



The integral over k in the previous equation then yields 



1 



a b2 = ^ooe I dx I dy 



while the integration over the variable y gives 



1 



v/(l - x)(l -x-y) (Vx(l -x) + y/yj 



a b2 = ^o'oe / dx 



7T + 2i 



1 -x 



■ sec \/x 



Finally, after performing the remaining integral over x, we obtain 



4 - 7T 



Cfc2 



a e 2 . 



(G18) 



(G19) 



(G20) 



(G21) 



Similarly, using Eq. (|G17|) after shifting p — xk > p, and integrating over p, the term crj^ in Eq. (|G12|) becomes 



-(2) 

a bl 



-^a e 2 n 3 f dx f dy(l-x- V y 1/2 x 
16tt Jo Jo 



1/2^-1/2 



d 2 k 



x(l - x)k 2 + y (§)' k 2 + (%Y 



(G22) 



Integration over the remaining momentum variable yields 

el rl-x 



~(2) 1 ill 

(7^=-— <r e I dx 



dy 



^x{\ -x-y) (y + yjxy(l - x)J 



~ (2) 

After integrating out the Feynman parameter y, the term d bl is 



, ( nl , sec 1 v r 



~(2) * 



7T 2 

= --eroe 



Let us now calculate the term <7b3 in Eq. (|G9[) . Using Feynman parametrization (|E16I) . we can write 



(1 - x - y) l l 2 x- 1 ' 2 



(p - xk) 2 + x{l - x)k 2 + y 



3 • 



(G23) 



(G24) 



(G25) 



After shifting p xk — > p, and integrating over p, we have 



±-a e 2 nj dxj dy {l-x-yf^x 3 ' 2 
d 2 k fc(ft 2 + 8fc 2 ) 



(27T) 2 



,{l- x )k 2 + y(%) 2 fc 2 + (§) 



while integration over k then gives 

1 



Cfc3 = —aoe 

S7T 



^1-x-y 


x(l - x) + 4^xy(l - x) + 2y 




{l-xf/ 2 (y+ y /xy{l-x) j 





(G26) 



(G27) 



which after calculating the integral over the variable y becomes 



1 



a b3 = — cr e 2 / dx (1 - x) 3/2 2(x - l)\fx + 7r(l - x 2 ) + 2i^x 3 (l - x) sec 1 yfx 
8vr J L 



14- 3tt 



24 



-cr e 
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(G28) 



Therefore, we calculated all the terms needed to obtain 
the Coulomb correction to the conductivity within this 
regularization scheme. Using Eqs. (|G16|) . (IG21|1 . (|G24I) . 
and (|G28p , we obtain the final result as found in Ref. |9| 



v {c) = a a +a b 



, ^( 2 ) 



bl ~^°bl ■ ■ -™ ^ 

(G29) 

which thus yields C = (25 — 67r)/12, different than one in 
Eq. (|82[) obtained using dimensional regularization with 
Pauli matrices in D = 2 — e spatial dimensions. 



25 - 6vr 



-cr e , 



spatial dimensions, provided extra care is taken in eval- 
uations of the integral that defines it. The result is then 
in agreement with the general dimensional-regularization 
scheme used in the rest of the paper. 

As suggested by MishchenkoiS, the expression for the 
first-order correction to conductivity may also be conve- 
niently written as a sum of three terms 



a = cr 



(HI) 



where 



Appendix H: Direct evaluation of the a. c. 
conductivity from the Kubo formula in two 
dimensions 



In this section we show that the value of the coefficient 
C may also be unambiguously computed directly in two 



e 2 w 



fpfk 



Up - k cos 6> pk 



2 , d^k 4 

°b = e ^ I Vp-kCOSflpk 



lo 2 - Ap 2 
p 2 (u> 2 + 4p 2 ) 2 

(uj 2 /pk) cos 9 
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{lj 2 +4fc 2 )(w 2 +4p 2 )' 
(H3) 



-2e 2 w 



d 2 pd 2 k k—p cos 9 

Vp-k cos pk 2 . -. 2 - 2V (H4) 



(2tt)4 



p 2 k(uj 2 + 4p 2 ) ' 



Here we have set the fermi velocity to unity for simplicity, 
and taken uj to be the Matsubara frequency. 

Now we show that all three terms are UV convergent, 
and when summed yield the same result as the dimen- 
sional regularization. First, 10 



<Jn + (r' a = 2e 2 co 



where 



d 2 p 



(2tt) 2 p(l + 4v 2 p 2 ) 



v F 1 



e 2 A 
— In — 

4 p 



0(V 2 ), (H5) 



(H6) 



is the renormalized velocity. A simple change of variables 
then gives 



where 



Me 2 ) 



dv„ 



v p dln(p) 



+ 0(e 4 ) 



(H7) 



(H8) 



is the beta-function for the velocity. The coefficients in 
the series expansion of (3 v (e 2 ) are universal numbers, and 
therefore, 



O"0 



(H9) 



in agreement with ref. The result also agrees with the 
brute force numerical integration, in which the integral 
over the angle is computed first and exactly, to be fol- 
lowed by the UV-convergent integrals over the momenta, 
which are then computed up to a large cutoff. 

The second integral is completely convergent, and was 
solved by Mishchenkc-^, with the result 



(TO 



e 2 . 



(H10) 



We have also reproduced this value numerically within a 
tenth of a percent. 

Finally, the third integral may be written as 



/ 



2 e 2 

— 7T hm 



Ai/u 



dp 



p(l + 4p 2 



A 2 /u 



2n 



d6 cos 6(p 2 + k 2 -2pk cos 0) 



1/2 



(Hll) 



r 



where we have carefully retained finite upper cutoffs on 
the momentum integrals. One finds 



2e 2 



Ai/u 



dp 
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— = 7t lim / — : — „, 

cr 7T 2 u-X)J p(l+4p 2 ) 

where 

2\x-y 



F p, 
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F(x,y) = 



3xy 
{x 2 + y 2 )E 



{x + yfKl- 



Axy 
(x - y) 2 J 



Vh13) 



Axy 



(x- y y 



and K(z) and E(z) are the elliptic functions. The sin- 
gle remaining numerical integration quickly converges to 
a value quite independent of the ratio A1/A2, and we 
obtain 



— = 0.2498e^ -> — , 
cr 4 



(H14) 



where the last equality is the conjectured exact result. 
All put together gives 



a 
co 



4 7T 

3^2 



If 7T 



e 2 , (H15) 



in agreement with the procedure of dimensional regular- 
ization. 

It is instructive to see how in this calculation the value 
C = (25 — 67r)/12 would arise. If we follow the usually 
safe practice and take the UV cutoffs to infinity before 
all the integrals have been performed, in the present case 
the result turns out to depend on the order of integration. 
In this way performing exactly the integral over k first 
we find 



/ 

(TO 



2f_ 

' 7T 2 



dp 



p(l + 4p 2 ) 



{-up) 



(H16) 



which ultimately yields C = (25 — 67r)/12. Computing 
first numerically the integral over p and then the remain- 
ing integral over fc, on the other hand, leads to e 2 /4 in- 
stead. Of course, by the very definition of the integral, 
the correct way is to take any limits of the integration 
bounds only after all the integrals have been already per- 
formed. It is pleasing to see that this then leads to the 
same result that is obtained by the general dimensional 
regularization, which is constructed so to preserve the 
crucial symmetries of the theory. 
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